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ABSTRACT 

U»Uke the signal received over a line-of-sight path, scatter signals arrive at the 
receiving site from a continuum of directions »ith iutens|tj^es that may be 
by a directional pattern similar to an antenna pattern. This report shows how the 
mean signal power available at the terminals of a rece.vmg antenna may be ex- 
pressed in terms of the antenna pattern and of a pattern of mcomiag power density. 
Lthods of measuring the power-density pattern at the rece.vmg site are 
Designs maximizing the available power are obtained for two rather general types 
of antenna when they are illuminated by arbitrary power-dens.ty patterns. Numerical 
results, obtained for an assumed Gaussian-shaped power-seos.ty pattern, suggest 
that only a very small increase in available power may be obtained by readjusting 
the pattern of an antenna that was initially adjusted to have maximum plane-wave 
gain in the direction of maximum power density. 

A large array of small antennas may be subject to “gain loss” but, if the 
tions between array elements are allowed to vary in time at the fading rate, the 
array will not exhibit “gain loss." Appropriate electronic circuitry for the inter- 
connections is discussed. 
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TOPICS IN THE DESIGN OF ANTENNAS FOR SCATTER 


I. INTRODUCTION 


Practical communication links utilizing the weak, fluctuating scatter signals received from 
a remote transmitter have been demonstrated.^ Unlike signals received over a short line-of- 
sight path, these scatter signals arrive at the receiving site from a continuum of directions with 
intensities that may be described by a directional pattern similar to an antenna pattern. This 
spread in the directions of arrival of the signal components gives rise to a phenomenon that has 
misleadingly been called "gain loss"; The strengths of scatter signals from a remote trans- 
mitter, received on adjacent, dissimilar antennas, are not in the same proportion as the plane- 
wave gains of the antennas. 

As suggested by Schott,^ the mean signal power available at the terminals of a receiving 
antenna is proportional to the sum (or integral) over all directions of the pattern of incoming 
power density weighted by the antenna pattern. In arriving at this result, Schott makes the 
plausible assumption that the plane-wave components of the signal arriving from, different direc- 
tions are uncorrelated. This assumption, which is also made here, is implicit in the concept of 
a pattern of power density; the assumption is discussed in greater detail in the next section. 

The aim of this report is to provide a method for arriving at a "best" antenna design for 
the reception of scatter signals. Certainly the "best" antenna should deliver the "greatest" 
signal strength to its associated receiver but, since scatter signals are subject to fading, the 
choice of a "greatest" signal implies a criterion of measurement. Also, in general, the larger 
the antenna, the greater the signal it will deliver. In the sense of this report, the optimum re- 
ceiving antenna is that stationary structure of a given size which delivers maximum mean signal 
power to a matched load. 

Although the pattern of incoming power density gives necessary and sufficient signal data 
I for the design of an optimum receiving antenna, the shape of this pattern depends on the shape 

of the transmitting antenna pattern as well as on the behavior of the scattering mechanism. In 
other words, a detailed description of the scattering mechanism itself would be necessary for 
' the simultaneous optimization of both transmitting and receiving antennas. Such a description 

I is beyond the scope of this report. Nevertheless, the remainder of this introductory section is 

j given to an elementary discussion of the effects that the transmitting and receiving antenna pat- 

I terns may have on the received signal power. ^ 

Following a simplified version of the Booker -Gordon model of the scattering mechanism, 
let us suppose that the transmitting cintenna emits power uniformly in a circular cone having a 
small angle 0^, as in Fig. 1. Similarly, let the receiving antenna accept power uniformly only 
in a cone having a small angle 0^. Suppose that the region above the surface of the earth is 
entirely uniform except for some subregion having a volume V, where scattering takes place. 
Suppose that each volume element, AV, of V scatters power isotropically, so that KaV watts 
per unit solid angle are scattered when the volume element is illuminated with unit power density 
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from the transmitting antenna, and let K be very small so that V intercepts negligible power, 
and multiple scattering may be neglected. Let the linear dimensions of V be small compared 
with the distances from V to the transmitting and receiving sites. Finally, let V be entirely 
contained within the common volume defined by the intersection of the transmitting- and 


receiving-antenna cones. 




Fig. 1 . The scattering region is contained in the 
common volume. 


Fig. 2. If beamwidths are reduced, the common 
volume is contained in the scattering region. 


The power gain of the transmitting antenna is inversely proportional to 9^ , and the power 
gain of the receiving antenna is inversely proportional to 9^. The signal power available at 
the terminals of the receiving antenna is proportional to the product of the antenna gains, the 
scattering coefficient K, and the volume of the region in which scattering takes place. V. As 
far as the antennas are concerned. 


(received power) 


If the angular width of the receiving antenna pattern, 9^, were varied, the received power would 
be just proportional to the gain of that antenna: In the present elementary circumstance, the 
"gain loss" phenomenon does not appear provided the scattering region V lies entirely within the 
common volume. In other words, the dependence of received power on antenna gains is the same 
as that for a line-of-sight path if the antennas cannot "resolve" the scattering region. 

On the other hand, if the antenna beamwidths are reduced to increase the received power, 
a point is reached where the common volume no longer contains the scattering region.V. Let 
us next examine the dependence of received power on antenna beamwidths when this transition 
has been completed; that is, let us assume that the common volume of the antenna cones is 
entirely contained within the scattering region V, as in Fig. 2. Suppose that the transmitting- 
antenna cone is the larger, and that it is pierced through its center by the receiving-antenna 


The region in which scattering takes place is now just the common volume, which can be 
described roughly as a cylinder having a height proportional to 9^ and a cross-sectional area 
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of this cylinder is proportional to 9^ 9^^ the received power 


proportional to 9^ . The volume 
is proportional to it and to the product of the antenna gains; 

(received power) 9^ 9^ ^ 

No further increase in received power can be obtained by reducing the receiving-antenna beam- 
width 9 , and the received power increases at only half the rate of the previous case when the 
transmitting-antenna beamwidth 9^ is reduced. Both antennas seem to exhibit gam loss. 

Suppose that the transmitting-antenna beamwidth is reduced to increase the receive power. 
When the two beamwidths are about the same size, the received power will cease to increase 
because the assumption about the shape of the common volume will no longer be 
the transmitting-antenna beamwidth 9^ is the smaller, the common volume is more nearly p 

portional to 0^, 

(received power) oc - — 

Bu. .. .his poist. receiv.ng -a...h.. hs.n..l» n».y h. --e... th. r.celv.d 

power. In en, case. 1. appears *•., I" g«ne™l. the most eeonomle.l design ot a pr.C.ea 
scatter link would call for equal transmitting- and receiving-antenna beamwidths. 

Before closing this discussion of the joint effects of the transmitting- and receiving-antenna 
patterns on the received signal power, it is appropriate to discuss the information required for 
L simultaneous optimization of both antennas and to remark on the difficulties of obtaining this 
information by radio measurements. In general, each volume element m the scattering region 
can be expect^ to scatter power according to some pattern, generally not isotropic, and the 
shape of this pattern may depend on the position of the volume, ^lilement in the scattering region, 
we may regard the magnitude of the scattering pattern in the direction of the receivmg an enna, 
produced when one volume element is illuminated by unit power density from the tr.nsr..m g 
antenna, as the scattering efficiency of that volume element. Thus, even if questions of polari- 
zation are somehow set aside, it would be necessary to know the scattering efficiency throughout 
the scattering region in order to calculate the received power when given transmitting an re- 

ceivine antennas are used. , 

A measurement of the power density arriving at the receiving site gives only two- 

dimensional information about the scattering region; the power density is a 'inc ion o e 
azimuth and elevation angles. The pattern of arriving power density is like a photograph of 
the illuminated scattering region as seen from the receiving site and is not enough irformat o 
to determine the behavior of each point in the volume. The additional ^ 

be made available by a power-density measurement at a second receiving site does ^ot appear 
to be helpful because the scattering pattern of each volume element would contribute differently 

to the new "photograph." 

The remainder of this report will deal with the density of power arriving at a single re- 
ceiving Site when a fixed transmitting antenna is used. As shown in the foregoing discussion, 
the results will pertain only to the effectiveness of the receiving antenna in extracting power 
when it is illuminated by this power-density pattern. 
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II. CROSS-CORRELATION BETWEEN ARRIVING SIGNAL COMPONENTS 

In this section, we will consider the assumption implied by the introductory discussion, 
namely, that the plane-wave components of the signal arriving from different directions are 
uncorrelated. At the receiving site, the electric field due to the arriving signal rnay be resolved 
into an angular spectrum of incoming plane waves. The antenna weights the amplitude of aplane- 
wave component F. arriving from a particular direction i, in accordance with the amplitude of 
the antenna pattern in that direction, a., and the open-circuit signal voltage V at the antenna 
terminals is the vector sum (or integral) of the weighted plane-wave components: 

V= 2-1^1 • 

The index i is associated with the direction of arrival of the i^'’ component. , ^ . 

The amplitudes and phases of the plane-wave components F., and hence the amplitude and 
phase of their weighted sum V, vary in time as a result of motions in the scattering region and 
possibly also because the transmitted signal is modulated. Temporarily, it will 
to assume that the transmitter is not modulated, so that attention is focused on random luctu - 
tions caused by the scattering mechanism. Later we will show that the results also apply to 
modulated signals. .r.»=n<;that 

I„ mi. report, the not.tio,,. Av.( ) »d O r.pre.... a Ome average, aad )* 
me comple. eoaiag.te 1. » be tmten, U.lng tbem, m.e m,.n-.,uare terminal ' 

proportional to the time-average available algn.l power at me r.oe.v.ng-antenn. termmal., may 
be written as » 

VF-Av.{(ya,Fo^y.gFg)‘}-2 2.iV-hF; . 

The cross-correlations cannot be r egarded as power densities because although their 

weighted sum over ml direcUons (i and k) is proportional to the total power, it is not possible 
to associate a single direction of arrival with F. F^*. (Is the direction i or k?) On the other 
hand, if 

F. F * = 0 for k ^ i 

the mean-square terminal voltage could be written 

and a single direction of arrival i could be ascribed to each term contributing to the total power. 

But the above assumption places too severe a restriction on the arriving signal components. 
It is neither reasonable nor necessary to assume that two signal components, arriving from 
directions that are arbitrarily close together, are uncorrelated. A physical antenna, havmg 
finite aperture dimensions, has a smooth pattern that does not change abruptly with small 
changes in direction: for radian changes in direction less than the reciprocal of the largest 
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aperture dimension expressed in wavelengths, the antenna pattern remains essentially constant. 
Thus, depending on the aperture dimensions, there is a "cone" of directions within which 
a « a. . Let it be assumed that 

k 1 . . 

F. = 0 for all k such that 3.^4 


The mean-square terminal voltage can now be written 

i k 

and the second sum may be regarded as the power density, since a single direction of arrival 
i may be ascribed to it. 

The relation, (2), and the resulting simplification, (3), indicate the nature of the assump- 
tion of incoherence that was implied by the introductory discussion and that will be made in this 
report. Essentially, we assume that no antenna we would consider building will have a narrow 
enough beamwidth to enable the antenna to resolve the correlation angle: the increment in angle- 
of-arrival over which arriving plane-wave components are appreciably correlated. Note that, 
although we assume no correlation over angular spreads as large as or larger than the beam- 
width of some "biggest" antenna, we make no assumption about the shape of the c 
for angular increments less than the correlation angle; none is necessary. It is 
able, for example, that the cross-correlation would decrease more rapidly with 
elevation than with increments in azimuth. 

Since no experimental evidence tending to invalidate the assumption of incoherence made 
here has come to the attention of the author, the assumption will be considered in the light of 
various theories of the scattering mechanism. In terms of the Booker-Gordon theory, we 
assume that the angle at the receiving site subtended by a length in the scattering region equal 
to the scale of turbulence is less than the beamwidth of the "biggest" antenna. Airborne refrac. 
tometer measurements, as interpreted by Gordon,'^ yield values of the scale of turbulence of 
the order of 50 meters for tropospheric scatter. Even for a short 100-km path, for which the 
distance from the scattering region to the receiving site is about 50 km, the correlation angle 
is of the order of 


r OSS- cor relation 

quite conceiv- 


= 10"^ radian = 3.4 minutes of arc. 

50 km 

An antenna capable of resolving such a small correlation angle would have aperture dimensions 
at least of the order of 10^ wavelengths - a 30-meter dish receiving a 10,000-Mcps scatter sig- 
nal, for example. Thus, if the scattering mechanism is accurately described by the Booker- 
Gordon theory, the assumption of incoherence wUl be valid in most practical cases. Only in 
the case of large antennas, extremely high frequencies, and short path lengths, will the assump- 
tion become questionable. 

Bailey, et have applied the Booker-Gordon theory to E-layer scatter, using a value of 
100 meters for the scale of turbulence. Since ionospheric scatter becomes important at much 
lower frequencies (20 to 100 Mcps) and for greater path lengths (of the order of 10^ km) than 
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does tropospheric scatter, it is clear that the assumption of incoherence is in complete agree- 
ment with this theory of E -layer scatter. 

vm.™ .nd dis=«..in. .u.d»le„c= ,ro„ . „„,d-d,„»lc po... o ho. 

,ha. UioeUe energy rn the ..nroepher. i. dlesLpa.ed by being Iranelerred ” 

eddy to sdcceeslvel, .m.Uer eddl... In epUe ol .be respiting .s.or.men. of eddy sizes tb 
attribute the scattering of radio waves only to those eddies having a size o tie oi 
vertical wavelength, 

2 sin 0/2 

where A Is .he wavelength end 6 1. die sestterlng angle through whleh the radio ray. are l^n. 

r:;==r::— 

Zcy, but .e ape . pro 

,oL havl b. m .be order of S.O voters If the 
anthn. were to resolve hre correlation angle. Thus, In most practical situations, d,e assump- 
tion of incoherence is compatible with the Villars-Weisskopf theory. 

Perhaps the most interesting justification of the assumption of incoherenc - 
H .tration that in certain circumstances, the maximum signal power available at the ter- 

iTaTs of an antenna receiving coherent radiation is in^ P-portional to the antenna gam. 
minals an observed with properly pointed antennas. 

r:r."r=ri;'=t-= 

portional to ja/ and, if the indices are appropriately assigned to the directions, 

' ^ 2 
w = 2 la/ 


t onal to the total power radiated from the receiving antenna when it is excited with 
rs ^ proportion,, to di. radl.tlon rcs.st.noe K. ol the receiving 

since W .1/1 HI' R.. »d we h.ve chosen iihl. current 1. Thps, the signal power .v.ll- 
able at the terminals of the receiving antenna is proportional to 





from (II »d ( 41 . Explicit expression, corresponding to ( 4 ) mid ( 5 ) wlU he dewloped in Secs. HI 
and rV. 
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rn is confined to directions for which the arriving plane- 
dated; the sums in (5) will cover only those directions. 


"Wonf •xT/rr.r "II . p™po.«o... .« .e ^ 

We may write 


F.(t) = g(t) h with lg(t)l 


complex time function 
me functions. In that c 


; independent of directic 


mplex numbers f^^ 


and, by the Schwarz inequality, the numerator is less than or equal to 

E la^l" E Ifkl' ■ 

i k 

».». .1... «.P •• "'r- “ *’• " 

W, « X I'll' - X If^' ■ 

i i 

The InequaUty holds with the equal sign when a, = fi*. In that case, the more direction 
Usl 1 . die greater hi, antenna beamwidth, die greater will he die received power. If 
: rddZ; iSP . m die rimge o' direction, covered by die sum, die received po.erwotdd 

TalTZrs'Zt ItrZIuM "bl pCn.ed out .a. the antenna pattern adiu.tmen. 
a . f» required for thl. unusual conclusion will, in general, be practical only if the magnitude, 
and phase, of the fj do not change abruptly from one direction to the ne*t. 


SO far, we have iustlfied the assumption of incoherence, which allows the received s.gn^ 
power to be written in terms of die tuilenna pattern and of the p..«r. of Incoming power density, 
as in (31. Temporarily, «e assumed that the fluctuations of the plane-wav. components F, 
“Jly due to nuchiation. in the sc.t.ermg region. But if the shiq» of the power density pat- 
tern were affected by die additional fluctuation, caused by moduladon of the transmitted signal, 
die forthcoming mi.enn. design considerations, which are based on the power-density pattern, 
would also depend on the modulation. Fortunately, it appears that this is not the case, as we 
sZll now see Wide-hand scatter signals may be distorted because of multipalh transmission; 
it should be noted that we will not attempt to disprove the existence of this effect. 
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we wish to Show iha. die mcorning power density is proportional only to the mean transmuted 
no..* and 1. independent of the modulation. I. will be convenient to calculate die power receive 
power ana is inaei . , H n tvuical direction and to show that this power does not 

Witt, a very directive antenna, pointed typi „„qal.,lon. If this hypodi.tic.l 

SCATTERER p| receiving antenna is bigger than the "biggest'' 

E antenna we might consider building, and if it 
R has a beamwidth of the order of the correlation 
angle, the power density is proportional to the 
1 me an- square signal voltage at its terminals. 

surface then, that our hypothetical 

LV vm receiving antenna delivers an output only when 

TRANSM|TTmc RECEIVING illuminated by plane -wave components arriving 

from within a cone having an angle of the order 
Fig. 3. The received sign<il Isa superposition of signals correlation angle, and that the pattern 

scattered from the individual scatterers. ^^pUtude is constant within this cone. 

in what follows, it will be necessary to work with a model of the scattering 
We will assume that the signal voltage at the receiving-antenna terminals can e c^ "" J 
Te lum of signals scattered from moving scatterers in the receiving-antenna beam, as shown 
in Fig. 3. Neglecting the possible dispersion effects of a single scattering, the contribution 
the received signal voltage from the scatterer can be written 


v(t) = 2 bp i(t - y . 

«« “ • «; - — dei’;“:r:: 

the transmitting antenna. The time delay t changes as the 

u T dVvm ..tb qqpfqtterer and thence to the receiving antenna, t changes as me p 
mitting antenna to p „^fficient b is the ratio of the component of v(t) coming from the 

scatterer moves. The real coefficient b^ the _ 

scatterer to the transmitted signal measured seconds earlier, I I. ^ 

1 frth of the scatterer, the antenna patterns, and inverse distance, 
scattering strength of the p scatter signal as the projection of a rotating 

The usual expression for the instantaneous transmitted signal tne pr j 


.b„. w IS fbe usuim. e.rvl.r frequeuey. The vector ILl bescrlbe. 

Zorn process; If Umrsmltfer were umooqulat.q, Wl - cou..».. The received vol.sg. c„ 
be similarly represented. In term, of the vector eapr.ss.on for i(tl. 


hp"'>’x>>,.''"‘’'ic-v] 
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We will call the complex coefficient 



so that the time-dependent vector giving the magnitude and phase of the received voltage is 


V(t) = 2 C^ I(t - T^) . 

^ r does not depend on the modulation, 

According to the foregoing analysis, the magni scatterer again independent 

and its phase is the c^ phase delay of the path ^ , 1, ,,, ^.^Uude 

of the modulation. Actually, a single scattering may be on the signal fre- 

and phase of the scattering strength of a single scatterer ^ slowly with signal 

,uency. Also. f ^ 

frequency. These effects would make Idence^s for most communication signals, since 

we are certainly justified in neglecting these dependences f.equency. It 

they generally occupy a very narrow band of ^„tions in the 

r:;^rp„r4r -i. ~ — - v - 

“"w.p ..CV.. “ 

= 1/2 !v(t)l^ = 1/2 Av.fc C^ I(t - y E 

In V 


= 1/2 2 2 


and, in view of the statistical independence of the c's and the I's, 

7^) = 1/2 2 2 ^ • 

Finally, if the cross-correLion between the gains of paths through indi^dual scatterer.^ 
decreases essentially to zero for such short differences in time-of-fligh , ^ ,, 

the Llayed modulation is essentially unchanged, that is, if 
c = 0 for all v 

p V (fa) 

such that I(t - T^) ^ I(t - Tj^) . 
then the mean-square received voltage becomes 

7w../2W?2i:'v; • 

P V 

and the received power is seen to be proportional to the mean transmitted power and independent 
of the modulation, as required. 
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■ ■ t the last result it has again been necessary to assume a short correlation 
arriving at the • time-of-flight with distance: 

distance - in this case, in me u because of the usually small 

essentially in the vertical direction in ^eattering process, the time 

scattering angles 9 through which the radio ,,, ^th por 

of flight over the p^*^ path changes relatively slowly p,,^,.tly 

example. Villars and Weisskopf have between paths passing 

responsible for scattering) in sue a way wavelength. The difference in time-of- 

through the top and bottom of such an e y is j carrier. Thus, according to the 

flight over these paths is therefore equa to a per a bandwidth of the order of the 

Villars-Weisskopf theory, a scatter result. (7)'. 

carrier frequency to invalidate the (6) are generally 

The signal bandwidths required by e o j^icrowave frequencies and for the greater 

less than the carrier frequencies, particu ar y or question 

path lengths. Nevertheless, the bandwidths are great enough so 

about the validity of (6). available at the terminals 

.P.S ..con, «. Z Z, .L .... .... ...pr...... - 

of an antenna illuminated by scattere ra provided that plane-wave components of 

.nterpreted as the sum of weighted power and this assumption 

ibe si^al arriving from appreciably ^ent d^^^^^^^^^^ 

v... ...... ...» .. ....... - - «• 

this report. 
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HI. DESCRIPTION OF THE ANTENNA 


>plele descHpUons of the antenna pattern and of the pattern of incoming 

optimum antenna. The antenna is considered m this section, 

considered in Sec, IV. In arriving at explicit 


will be required for the design of 

s behavior when receiving scatter 


ignals 


r by the pattern of plane 



will be extended to three dimensions and used as fields it would produce if 

The receiving antenna will be characterized either by the eiecir 

of one amoere were applied to its terminals c 
a radio-frequency current of one ampere weie 

waves that would emerge from the antenna with the same 
excitation. Points in space will be located on a system 
of rectangular coordinates, with the antenna in the vicin- 
ity of the origin. As shown in Fig. 4, it will be conven- 
ient to speak of the positive z-direction as being "up," 
although the coordinate axes may always be oriented to 
suit a particular situation. It will also be convenient to 
assume that the region above the xy-plane is free space 
and that all the radiation from the antenna crosses the 
xy-plane and leaves the antenna in directions above the 
horizontal; at least, it will be assumed that any obstacles 
above the xy-plane are at great enough distances so that 
their presence does not materially affect the driving- 
point impedance of the antenna. Thus, we imagine an 
antenna located more or less at the origin of coordinates, 
radiating into the upper hemisphere, and having all its 
active conductors, supporting members, etc., tbe 
lower hemisphere. 

Although the restriction to one hemisphere - oimoler 

rrr.-; .« p-. --p 

separate solutions for each hemisphere. 

We will measure distances in wavelengths, so that 

j2ir(f t-z) 

Re{Ee } 

,h. .xpre..io„ for . coupon.,, of th. ol.ctric field of • pL". of o-_^. 

which hos l.K the .otenn. and 1. traveling toward the een.th. We .hall rep 

ponent as 

^ ^-jZTTZ 

it helhg under., ood'that , hi.’ and .Imll.r expr.e.lons for rwflo-fr.gu.noy quantitle. are to be 
multiplied by and the real part taken. 



Fig. 4. Geometry of the space surrounding 
the receiving antenna. 


Although di. r...rlc,lon to one hemlepher. hold. 
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In giving the direction of travel of a plane wave, we shall use the direction cosin 
n, which are the projections of a unit vector on the x-, y- and z-axes, respectively, 
Fig. 4. We regard the direction cosines t and m as the primary variables; we will ui 


;o write n in terms of I and m. Thus, i 
in the direction (f, m) as 


night represent a component of a plane \ 


-j27r<xl+ym+zv 1-f -m ) 


i mean that the square r 
o the upper hemisphere. 


o be taken with the posit: 


traveling into the upper hemisphere. However, we snaii - 
brevity. 

If the projection on the xy-plane of the intersection ol 
unit sphere is d« dm, the increment of solid angle itself is 


e sign, since tne wave is 
1 in place of the radical for 


-emental solid angle with the 


We shall often have occasion to integrate plane-wave components and other functions of direc- 
tion, (f,m), over all directions in the upper hemisphere. These integrals will be written as 


The second integral may also be taken to cover the inside of the unit circle, 1 + m - 1. 

Now we will use a physical argument to express the fields, produced when the antenna is 
excited with a current of one ampere, as a superposition of the plane waves that emerge from 
the antenna. Let the x-, y-: and z-cOmponents of the electric field at the point (x,y, z) be the 
complex functions, y, z) , Ey(x, y, z) and E^(x, y, z) . respectively. As shown in Fig. 4, let 

•Fjj(f ,m) dJl be the horizontally polarized component of the electric field of the plane wave that 
leaves the antenna and progresses into the incremental solid angle dSl in the general direction 
(l.m). The phase of Fjj«,m) will be referred to the origin of coordinates. Similarly, let 
F (f ,m) dSl be the vertically polarized component. As usual, we have taken the horizontal com- 
ponent to be tangent to a horizontal circle of constant latitude on the unit sphere; the "vertical" 
component is tangent to a longitude circle. 

At the point (x,y,.z), the x-components of the plane-wave components add to give 
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-j2T(xl+ym-tzn) 


Ey(x,y,z)= y 
E^(x,y,z)= ^ 


Using the index i to represent one of the subscripts x, y 
written as 




these three equations may all be 
.j27r(x{+ym+zn) ^ ^ (9) 


I 

i 


i 



The P's may be regarded as rectangular plane-wave components as contrasted with the F^s 

which are horizontal and vertical plane-wave components. It 

.L At- thp set of P’s differs in magmtude from the electric vecior 

Before continuing with the argument, we win pause 



n's are eliminated, there results 


(10a) 


-- X y '■ 

* ,1 + mean that the scalar product of a unit vector, pointing in the direc- 

electric vector is normal to the direction of propagation. This result is a consequence 
aving (correctly) taken the F's to be-normal to the direction of propagation ^ 

Suppose that we have two plane waves with rectangular components PilW.in) and Pi2(^ 

.d corresponding horizontal and vertical components F^.U. m) and F,,U, m)^ 
irst subscript in this case, k, on the F's to indicate the polarization. H or V. m the same way 
hat we have uied the first subscript on the P's to indicate which rectangular component, x, yor z. 
,s meant.) Using (iO). we find • 
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2 Pii Pi2 = 4 2 Fk; 


siPii' = 42 iFj 


terms of the rectangular plane-w; 
of (x,y) to a function of (i,m): 


n our attention to (9), which gives the field components in 
components. For a fixed value of z, (9) relates a function 




In the plane, z = constant, let E. have a Fourier transform, 
form also exists. The inverse transform could be written 


-jZTTzn ^-j27r(xf+ym) ^ 


But in the subregion, + have found that P! = R. It would appear that some "plane- 

wave components " which leave the antenna in the fictitious directions, i + m > 1, but which 
nevertheless may contribute to E., are missing from (9a). Indeed, the fields produced by a 
physical antenna on a plane in front of its aperture (z = constant) may vary so rapidly with posi- 
tion, (x,y), that the Fourier transform of these fields cannot possibly be zero everywhere out- 
sideV + m^ = 1. Let us examine the nature of the contribution to E. of these fictitious "plane- 


components." Writing nil 


^-jZTrixl+ym+z’Jt-P-m ) ^ 


Outside of + m^ = 1 the radical is imaginary, and we must reconsider its sign. If n were posi- 
tive imaginary, the contribution dEj(x,y,z) would increase exponentially with height z. which 
would not be compatible with the chosen source location. But if n is negative imaginary, the 
contribution decreases exponentially with height, which seems more^ reasonable. 

Since the contribution to the fields from the R (I, m) for which f + m > 1 becomes negligible 
at sufficiently great distances from the antenna, we associate the storage fields with these fic- 
titious "plane-wave components." The radiation fields, as we have already seen^ are ^accounted 
for by the components of those plane waves which travel in the real directions, 1 + m « 1. Thus, 
the complete expression for the fields is 
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and the direct transform, giving the rectangular plane-wave compone, 

H,,.^)=y]'E,(x.y,z)e^^’^^^^^-^^'^><ixdy . 


1 . m both of these expressions. Note that the fields at any point in the 
The height z IS a parameter in both of hese e p 

region above the xy-plane can be ca cu a e interrelated by (10a). The 

'-r:: 

::;:r rixu' „eu. ...» ......... ......ons wh„, 

IHl 1, .ppllc.bl,, apd .1... they match the boundary condl.lous .ppropr.a.o to .ha. rag.on, Th, 
Mauw... e,ua.l.ns appUcab.c to lb. .l.c.ric ,.cld In th. ,r..-ap.c. upper hem.spb.re are 
wave equation and the divergence relation, usually written 

£ = ep ^ . and V • ? = 0 . 

enre w. have m.a.ured diataueea In wavelength., the wav. equation become., <or the .Inu.oidal 


ax*^ ay sz 

•rvrto ar-f. a satisfactory approximation to the wave equation 

rh!rtLT,rrrrclT.“nrb.nL.ou. ampl..ud.. O, th. field 
luted wave. A. may be ..en by applying them to (.11, these expression. 

intend. The divergence relation may be written 


~:n" the upper heml.^.re, which ha. the xy-^an. as ... o.y houn^ry 

ar the source. The region will be excited by a physical rnitenn. located In the vicinity 
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origin „<! .omowha. b.lo. .he ay-plane. Surely, .he Fourier .rriiatorm of .he ti.l^ s P 
,h, x,-plan. by ,hie an.enn. ,111 eal.t, and the Inuer.e .ranef.rm will al.o eai.t. Thu. bound 
1'; condition, can be ntatched on .he xy-plane. The choice of .Ign for .he dlrec.ion coriue n 
given in (11), will msure that the fields given by (H) will decrease appropriately m intensity 
great distances from the origin. 


Next we will study the behavior of the fields at a point a great distance 


from the antenna and in the (real) direction m^, nj in the upper hemisphei 
point of space in question is 

(x,y, z) = (rf^, rm^,rn^) . 

We expect the fields at that point to behave like 


that the 


-jZxr 

Ej(rf^, rm^, rn^) ~ ^ — 


so it is appropriate to evaluate the limit. 




from (11). Since we are interested in the far fields, only those plane-wave components for which 
4 i will contribute, and the range of integration may be reduced to the upper hemisphere. 


r -327rr(/i +mm i-nn I 

L= lim re^^” \ nP.(i,m) e dS2 . 

r-*«> u h. 

It will be convenient to regard the exponent as the scalar product of two unit vectors: one, U^, 
having the direction {l^.Tn^.%). and the other, U. pointed at the incremental solid angle dfl: 


jZxrd-U-xfj,) 

L = lim r \ n P.(i,m) e dSl 

u.h. 

The geometry of this integral is shown in Fig. 5. We will integrate first around strips on the 
unit sphere in which the scalar product U • and hence the entire exponent, is constant; then 
we will add up the contributions of the strips. The strips are defined by circles of constant 
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latitude, formed with the polar axis in the direction of U^. 
the upper hemisphere, some of the strips are incomplete; 
defining the integrand to be zero below the xy-plane. 


Since the integral is taken only over 
we may account for this condition by 



Fig. 5. CT- lT is constant throughout < 
the strips. 



the 


Fig. 6. The strips are parallel 
Xp-plane in the new coordinate system 


The integration is facilitated by taking new axes, (X,p, r), in such a way that the directions 
of r and U coincide, as shown in Fig. 6. The new variables of integration will be the direction 
cosine r and the azimuthal angle 9 defined with the polar axis in the r-direction. We will define 

(n PU,m) 

Q(r 

With these changes, the limit lecomes 
L = lim 



Although we will not attempt an explicit expression for the function q(r), some of its proper- 
ties are of interest. The rectangular plane-wave components P,, which are the Fourier trans- 
forms of the components of the electric field set up on the xy-plane by a physical antenna, will 
be continuous, bounded functions of direction in the upper hemisphere. The function q(r) is an 
integral of one of the P. over a strip having a length that is also (in general) a continuous, bounded 
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function of thus, q(.) will be a continuou_s. bounded discontin- 

the particular case where the l^Zl ^ n-0. corresponding to the hor izon. 

|„„1 lUeiilL nP. decreases continuously to zero with 

,1 f t u at that point, making q(.^) again continuous_^and bounde .) 

I /■ ; Also, since we are interested in directions m the upper 

; ! hemisphere, the value q{-l) would involve only contribu- 

^^-^1 T tions from the lower hemisphere; in other words, q( 

-t A representative plot of |q(>') | is given in Fig. 7. 

Fig. 7. <\(y) is continuous and bounded. Consider the following integration by parts. 


J2r,r{i-v) ^ _ ej^’^^’qM) + i2;r: 

The second member on the right has been shown to be zero, 
to zero according to the Riemann-Lebesque lemma. Thus, 




.CO, the left-hand side goes 


0 = q(l) - 0 + jZTT lir 


- Q(l,e) de = jn„PiUo''"o* 
t Jn 


We have shown that 


which means that, at great distances r from the antenna, the electric field componeni 
direction (f , m , n ) are 


E.(rf^, rm^, rn^) SJ inoPi^o’^o' ' 


We note that the far field propagates in the direction away from the antenna, and that, since its 
components are proportional to the corresponding rectangular plane-wave components, the total 
electric field is normal to the direction of propagation. 

The far field may also be described in terms of the horizontal and vertical plane-wave com- 
ponents. Fjj and Fy. At the distant point (rl^, rm^, rn^), the horizontal component of the total 
electric field, which has the same direction as F^j in Fig. 4, is 
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a great distance, R = 7 


A being the wavelength, the Poynting v 


from the antenna. Its magnitude, 

s = 2 |eJ^ watts/meter^ 

(tj = characteristic impedance of free space) 

. power from fire an.enna oroa.lag . .,uare me.er of th. surface of a sphere " 

ms, centered a. the origin. The solid angle subtended by that se,u.r, meter at the origin 


K 

and the angular density of power from the antenna is 

^ - _S_ _ 2 |e. watts/solid angle 

-JIr 2 v , 1 

The concept of an angular density of power is useful for the purposes of this paper. When v 
refer to a power density, we will mean an angular density. 

using (13) and (10c), the density of power leaving the antenna and progressing in 

tion (1 , m) is 


I 
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and the total power emitted from the antenna is 
W = -^ 2 ^ |F^(f,m)|^ dfi 


= aI ^ „ |p.y,ni)|^ df dm 

‘ i u.h. 

= dmj! 


. r ,2 , 2 , 

since n is imaginary outside ot { + m = i. 

The ratio 

wU, m) 

W 

may be called the gain of the antenna in the direction (f, m): neglecting antenna losses, the ratio 
multiplied by 4jr is the gain relative to the gain of an "isotropic radiator." 

Before closing this section, we will digress to develop expressions for the transmitted power 
in terms of the near fields of the antenna. Starting with the last form of (ib) and using (12), 

w .|i SRejjjn <» dn, . 

Changing the order of integration, 

W =|i ^Relj]* E?(x,y,o) dx dy nP.(l,m) dm| , 

and the last integral is recogmzed from (11) to be 
. aE.(x;y, o) 


o that the transmitted power can be written 


2 f A rr 8E (x,y,o) 1 

- E Ji E?(x,y,o) dx dy^ 

i I .«> > 


which gives the transmitted power in terms of the fields and their z-derivatives on the xy-plane. 
The same result could have been obtained by integrating the real part of the vertical component 
of the Poynting vector over the xy-plane. Next we seek an expression for the transmitted power 
involving only the fields on the xy-plane. Starting again with (16) and using (12) to express both 
P. and P?, 
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W , ii jy „ d, dm E,(x,y,o| el"'«*>'”>dx d, 
i u.h. 

xy]’Ef(x',y',o)e-J'’^<"'‘+y'"’>dx'dy' , 


r with x' = X + u, y' = y + v. 


which becomes 


W= ndfdm^^ E .(x, y, o) dx dy 

i u. h. 

xfyEflxdu, ,*V. . 

X E.(x, y, o) E?'(x + u, y + V, 0) dx dy 


when the order of integration is changed. The last integral is the spatial auto-correlation of 
field components on the xy-plane: 


^..(u,v) =_^jE.(x,y.o) E:Mx + u, y + V, 0) dxdy - 

The second integral, which plays the role of a weighting function, is the Fourier transform of 
the real part of the direction cosine n: 
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Then let p = sin 0 : 


r ^ { I — 2 — 2 \ 2 

N^(u, v) = Ztt J^^2Wu +v sin<J)j sin* cos <f)d* 

2^^r(3/2) 



In terms of the auto- correlation of field components on the xy-plane, (18), and the weighting 
function, (19), the transmitted power may be written 


V = yy N^(u, v) <(). (u, v) du dv 
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IV. DESCRIPTION OF SCATTER 

It will be convenient to describe the scattered fields arriving at the receiving site in the 
same terms that were used in Sec. Ill to describe the fields that would be produced by the 
receiving antenna if it were excited with a unit current. The same coordinate system (Fig. 4) 
and similar notation will be used. The electric field components of the arriving radiation will 

be called E. (x,y, z). Similarly, the horizontal and vertical plane-wave components arriving 

from the direction (I, m) in the upper hemisphere will be designated m), and the corre- 

sponding rectangular plane-wave components will be P.^(f, m). The second subscript, a, will 
be used to indicate arriving quantities. The phases of the m) and of the P.^(f, m) (and 

also the amplitudes of these quantities when f ^ + m^ > 1) will again be referred to the origin 
of coordinates. 

We will assume that all the sources of radiation incident on the receiving antenna are 
located above the xy-plane. Physically, this radiation will be reflected (or scattered) from the 
receiving antenna and other obstacles below the xy-plane, so that the total field above the 
xy-plane will be a superposition of the incident fields and the fields reflected from below the 
xy-plane. The reflected fields will depend in part on the structure of the receiving antenna and 
its surroundings, but the incident fields will depend only on the transmitter and its antenna, the 
scattering mechanism and the geographical location of the receiving site. We will use only the 
incident fields (not the total fields) to describe the arriving scatter radiation. The field com- 
ponents E. , then, pertain only to the incident fields; they may, however, be regarded as the 
total field^components that would be present if the receiving antenna and other reflecting objects 
were removed from below the xy-plane, leaving free space. Similarly, the plane-wave com- 
ponents Fj^^ and P.^ represent plane waves arriving only from above the xy-plane; if they 
correspond to the total field, the lower hemisphere must be free space. This restriction to 
incident fields only will have to be considered again when the density of power arriving from 
the scattering region is to be measured in terms of the incident fields. 

The arguments of the last section leading up to and justifying (11) and (12) could be repeated 
almost word-for-word for the arriving quantities. The result would also be similar; the expres- 
sion for the arriving fields is 


Ei^(x,y,z)= yy P.^(f.m)ej2’^<^+y-+""Udm 


n is positive for i + m 1 

- .2 2 

n is negative imaginary for i ■ + m 
and the direct transform, giving the rectangular plane-wave components in terms of the fields. 


(21) 
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Note that the exponents in (11) and (21), and also the exponents in (12) and (22), have opposite 
signs. This is because the transmitted fields, which leave in the direction (f, m) are propagat- 
ing in the opposite direction from the arriving fields, which arrive frcp the direction (i,m). 

Malting the direction cosine n again negative imaginary for f + m >1 insures that the 
storage fields die out at the greater distances (more negative z this time) from the source. In 
fact, it will usually be reasonable to assume that the storage fields of the source do not contrib- 
ute to the fields arriving at the receiving site. In these cases, P^^(f, m) = 0 for I + m >1, 
and the range of integration in (21) may be reduced to the upper hemisphere. 

The relations (10) between the horizontal and vertical plane-wave components and the rec- 
tangular plane-wave components do not depend on the direction of propagation of tlie components. 
They are also valid for the arriving quantities. 

Up to this point, we have not mentioned the magnetic fields that must accompany the cor- 
responding electric fields. In general, since the magnetic field can be computed from the 
electric field with the aid of the appropriate Maxwell equation, it is not necessary to be con- 
cerned with the magnetic fields. However, magnetic quantities will be convenient in the forth- 
coming determination of the open-circuit voltage at the receiving-antenna terminals due to an 
arriving scatter signal. 

Temporarily, then, we introduce the arriving mjignetic-field components Hj^^(x, y, z) and 
their Fourier transforms Q.^(f, m), the (magnetic) rectangular plane-wave components. The 
magnetic field also obeys the wave equation and a divergence relation. By analogy with (21), 

(22) and (10a), 


Hia(x,y,z)= Q.^(i, m) di dm 


n is positive for I + m 1 
n is negative imaginary for t 


Q.^(i,m)= ff H.^(x,y,z)e-i2-(W+zn)^^,y 


i + m Qy + n = 0 


(24) 

(25) 


The similarity here is so complete that no comment is required. 

Next, we relate the arriving magnetic fields to the arriving electric fields, using the Max- 
well curl equation which is usually written 

VX?=-^ • 

In the present circumstances (sinusoidal steady state, free space, and distances measured in 
wavelengths), this equation may be written 
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- - ^ = -^^TTT) H^. 


Z!L _ = -i2i:7) H , , 

9x ay ■' za 

n terms of the arriving rectangular plane-’ 

m P_ - n P.„ = -7] 1 


e components, from (21) and (23), 


Pxa - ^za = V 
i P - mP = -7? 


It is enlightening to summarize in vector notation. Let P^(i, m) be an electric vector having 
the components P,,(i,m), P^^(i, m) and P^^(i. m); let Q^(i,m) be the corresponding magnetic 
vector; and let U(i,m) be a'vector in the direction (i,m). The wave equation requires that 

U • U = 1 , 

so that U is a unit vector. The divergence relations may be written 


U • P = 0 




U • Q = 0 , 

and the curl equations become 
U XP = -7) Q 


(10a') 

(25') 

(26') 




7j U X Q = P 


Although these vector equations have meaning only for the real directions, their algebraic 
counterparts, written in terms of the direction cosines, are also valid for the fictitious direc- 
tions f 2 + m^ > 1, as we have seen. In view of the definition of the Poynting vector, it is clear 
from either of the curl equations that PJf , m) and Q^(i, m) represent radiation that is propagat- 
ing in the opposite direction from that of U. In any case, the signs of the vector products would 
have been reversed if we had been describing the transmitted fields. It may also be seen that a 
number of these vector relations are redundant. 
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Now we are ready to develop an expression for the open-circuit voltage V at the receiving 
antenna terminals due to an arriving scatter signal. As suggested by Fig. 8, we will calculate 
the open-circuit voltage at the terminals of a small test dipole when the receiving antenna is 



Fig. 8. The current source and voltmeter may be interchanged 
without affecting the voltmeter reading. 


excited with a radio-frequency current of one ampere; then we will use reciprocity to infer 
that the same open-circuit voltage would appear at the terminals of the receiving antenna if, 
instead, the test dipole carried unit current. Let the test dipole be located at the point (x,y, z) 
above the xy-plane, let it be oriented in the i-(x-, y- or z-) direction, and let it have an effec- 
tive length of a meters, short compared with the wavelength. The open-circuit voltage at its 
terminals, caused by excitation of the receiving antenna, is then 


V = aE.(x,y,z) , 

and the same voltage would appear at the receiving- antenna terminals if the test dipole v 


excited with one ampere. 

On the other hand, the scatter fields arrive from a variety of directions, producing very 
complex field patterns at the receiving site. A single test dipole could hardly be expected to 
reproduce such patterns. Tentatively, we suppose that a current sheet on some plane, z = 
constant, above the xy-plane is capable of producing field patterns that duplicate the scatter- 
field patterns in the vicinity of the antenna. The current sheet may be considered to arise from 
an infinite array of incremental test dipoles, oriented in the x- and y-directions on the plane 
z = constant and excited with different currents. Since we know the open-circuit voltage at the 
receiving-antenna terminals due to the excitation of each of these dipoles, these contributions 
may be superposed to obtain the open-circuit voltage due to the current sheet. 

Let the components of the linear current density on the sheet be J^(x, y, z) and Jy(x,y, z) 
amperes/meter. Consider an incremental rectangle bounded by {x,x + dx) and (y,y + dy) in 
the plane z = constant. The total current flowing in the x-direction in this rectangle is 
(Ady) J (x, y, z) amperes, and the effective length of this rectangle of current is Adx meters. 

Its contribution to the open-circuit voltage is therefore 
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dV = - (Adx) E^(x, y, z) (Ady) J^(x, y, z) 


(The negative sign arises from the polarity conventions shown in Fig. 8.) Superposing the con- 
tributions from all the x-components and all the y-components of current density on the siee , 
we obtain the open-circuit voltage due to the current sheet: 




The Bio, -Sever. 1.. .he. .he megne.omo.lve lore, .round . deed loop equals ttie 

net current enclosed; directions are given by the right-hand rule. If the loop lies m the 
yz-plane with horizontal sides just above and just below the current sheet, as shown in Fig. 9 
L current involved will be related to J^. If the loop has differential length dy and vanishing 
height dz « dy, so that displacement currents can lie neglected, we have 
(Ady) J^(x, y, z) = (Ady) [Hy (x, y, z) - H +(x, y.z)] . 


= H ■ - 


where as in Fig 9, H+ and H ' are the y-components of the magnetic field just above and just 
LTow the current she^t. In Lder to compare the fields due to the current sheet with the 
incident scatter fields. 



imagine that the re- 
ceiving antenna and other reflecting objects are 
removed from below the xy-plane, leaving free 
space. In this case, the tangential magnetic 
fields above and below the current sheet are 
equal in amplitude and opposite in direction: 


o that J = 2H 


The corresponding relation. 


, . for the other components could be obtained by 

Si .PP.,..g Bio.-Sev„. ..w .o . .oop 


tangential t 
tenna is 


V = 2A^ (Ey H; - E^ H -) dx dy 


At this point in the discussion, the region below the current sheet is free space. Accord- 
ingly, H- and H- must be components of a magnetic field that is propagating mto the lower 
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-plane. In terms of the magnetic field 

plane just below the current sheet, the open-circuit voltage of the receiving an- 


V 


UNCLASSIFIED 



hemisphere. The corresponding plane-wave components are given by (24), and the missing 
vertical component is then given by (25). The magnetic fields anywhere below the current sheet 
could be calculated from (23), and the corresponding electric fields could be calculated from 
(26), (10a) and (21). In short, the fields anywhere below the current sheet are completely 
determined by the tangential magnetic field components and on a plane just below 
the current sheet. Since this tangential magnetic field is proportional, in turn, to the density 
of current on the sheet, we conclude that the current sheet is capable of reproducing the field 
patterns of the arriving scatter radiation at the receiving site, and we deduce that 


V = 2A^ (Ey - E^ Hy^) dx dy (27) 

is a general expression for the open-circuit voltage at the terminals of the receiving antenna in 
terms of the arriving magnetic field and the electric field produced by the antenna when it is 
excited with one ampere. 

Using (23) to express the magnetic field components in (27), we have 


V 




<^y “^xa - ^x V df dm , 


and, after reversing the order of integration, we recognize the second integrals as (12). Ac- 
cordingly, 


V = 2a2 


§ 




Using (26) to write the magnetic plane-wave components in terms of the electric plane -wave 
components, and using (10a), the open-circuit voltage can be written more compactly in terms 
of electric plane -wave components only; 

V = ^ 2 yy n Pj(i,m) Pi3(i,m) dl dm . (28) 

A word of caution about (28) is in order. This relation is valid even if the sources of the 
arriving radiation are near the receiving antenna; it is valid for the current sheet, for example. 
But when the current in the current sheet is turned off so that the pattern of the receiving an- 
tenna can be measured, nothing but free space is left where the sheet had been; the current 
sheet does not disturb the field pattern of the receiving antenna. A physical source near the 
receiving antenna, on the other hand, would generally reflect some of the fields from the an- 
tenna, which would invalidate (28). In that case, (27) would be valid if E^ and E^ were inter- 
preted as components of the total field on the plane z = constant when the receiving antenna is 
excited with unit current and the nearby source is turned off (voltage sources short-circuited 
and current sources open-circuited). 


Lm 
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However, we shall assume that the sources of the arriving scatte 
iistances so that they do not affect the fields produced when the recer 
shall also assume that the storage-field region of these sources does 
site. The scatter fields, then, will arrive only from the real directic 
integration in (28) may be reduced to cover only the upper hemisphere 

V = ^ n P.(f, m) P.^(f, m) di dm 

= ^ Fj^(i,m) Fj^^U,m) dO 


Having developed expressions for the open-circuit voltage at the receiving-antenna termmals 
due to an arriving scatter signal, (29), and for the total power radiated when the antenna is driven 
with unit current, (16), we can now write an expression for the time-average signal power avail- 
able at the antenna terminals. In terms of the unit current and the radiation resistance of the 
antenna, the radiated power is W = l/2 |ll^ R^, and, neglecting antenna losses, the mean avail- 
able signal power is 


• |j 2 y F.(f,m) F.^(f,m) dJlj I 

2 y lF.(f,m)l2 dJl 


Next we manipulate the numerator N of this expression for mean power, making use of the as- 
sumption of incoherence discussed in Sec. II. It should be noted that only the arriving quantities 
Fjj and F.^.^ are time-dependent; the time average will apply to them. Using a second, similar 
set of directional variables such that 


^=2 2 ^ F.(f,m)d0 j F*(f'.m') F.Jf,m)F^yf',m')dJ2' , 
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N = E E y F,(i, m) dn 
i k u.h. 

X F * (i + X, m + fi) F.^U. m) F^^(i + X, m + |i) . 

The second (double) integral involves the cross-correlation between arriving plane-wave com- 
ponents. The angular separation between the directions of arrival of the two components is 
given by tlie variables of integration X and p, which are the differences in the direction cosines. 
As in Sec. II, we assume that the cross-correlation 

F.aCi.m) F*^(i + X, m + p) 

will be non-zero only for such a small range of directions of arrival (X, p) that the antenna pat- 
tern will be constant over that range: 

F*(f + X, m + p) =!F*(f,m) . 

The direction cosine n', which generally varies more slowly with direction than the antenna 
pattern, will also remain constant over the range of appreciable contribution to the second 
integral. Finalty, the range of integration, which strictly speaking is a displaced unit circle, 
may be extended without changing the value of the integral. With these changes, the numerator 
becomes 

N = E 2 y F.{f, m) F^(i, m) dfl 
i k u.h. 

^ n F.gd, m) F^^d + X, m + p) dX dp , 
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Suppose that the antenr 
the available signal power is 


.lly-polarized waves when i 


..h. 

y lFjj(i,m)|^ df 


Let the power -density patter 


change to any other pattern of the same "area" would re- 
v ) duce the numerator integral, leaving the value of the de- 

nominator integral unchanged. Of course, such an opti- 
mization would require unreasonably large antenna- 
aperture dimensions. We will discover that, in general. 
Fig. 10. The available power is maximized if the aperture is specified, the optimum antenna has 
by the antenna having a pattern that is almost, .but not quite, maximum plane -wave gain; its 
sharply peaked in thedirection of maximum 

power density. pattern will be clustered about the direction of maximum 

power density. 

The case of antenna noise is also of interest. If "hot spots" in the sky are neglected, an- 
tenna noise may be considered to arise from a uniform distribution over the upper hemisphere 
of incoherent sources of all polarizations, such that the power densities are 
$ (i,m) = $ „(i,m) = , a constant 


The available noise power 


which is independent of the antetma pattern . Apparently, the anterma pattern that maximizes 
the mean signal power will also maximize the signal-to -noise ratio, since the noise power will 
not be affected by the pattern adjustment. 

We have finally arrived at the core of the problem. From the point of view of the antenna de- 
signer, the power densities are necessary and sufficient data for the design of an optimum antenna. 
His job is to produce an anteima having a pattern, determined bytheF.(f, m)t which maximizes the 
available power, (30), subject to appropriate restrictions. The proprietor of the receiving site, 
on the otherhand, will describe the power densities to the designer and specify some of the restric- 
tions. Section V will be given to a discussion of the power densities and suggestions for their meas- 
urement; the following sections will deal with some of the antenna designer's problems. 
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V. DETERMINATION OF POWER DENSITY 


It should be emphasized that tom power densities, given by (31) in terms of the arriving 
plane-wave components, were required to express the mean signal power available at the 
receiving-antenna terminals, (30). Since the numerator of (30), being a squared magnitude, 
is a real, non-negative number for any antenna pattern (for any set of F.'s), it must be that 


that the self-densities $ 
power densities 


and * are real non-negative quantities, and that the matrix of 


$ U,m) 


is Hermitian, with a non-negative determinant. For a particular direction, U.m). the deter- 
minant becomes zero if and only if the polarization of the plane wave arriving from (i.m) re- 
mains fixed in time as the signal fades. In this special case, the incoming radiation could be 
described in terms of two power densities but, in general, three are required; it wiU be con- 
venient to retain all four to preserve symmetry. If the available power had been written in 
terms of the rectangular plane-wave components, nine rectangular power densities would have 
resulted; these could, however, have been reduced to the above four with the aid of the rela- 
tions (10). 

It is the purpose of this section, then, to suggest means for determinmg the four power- 
density patterns: the four average quantities that describe the arriving scatter radiation. At 
the outset, it should be noted that, since the recorded strengths of received scatter signals 
generaUy show seasonal as well as diurnal trends, certainly the amplitudes and very possibly 
also the shapes of the power-density patterns may show corresponding long-term variaUons. 

If the patterns are to be used in the design of an optimum fixed receiving antenna that is ex- 
pected to deliver maximum mean signal power for its useful life, the power densities should be 
averages that are taken over periods of time long compared with the periods of variation of the 
received signal strength. Averaging periods at least of the order of several years are indicated. 

On the other hand, the proprietor of the receiving site might want mstead an antenna that 
delivers maximum mean power during periods of weak received signal. In that case, he shou d 
average the power densities only during such periods. Of course, the power-density patterns 
characteristic of weak-signal periods might show more statistical regularity, so that shorter 
averaging times could be employed. Nevertheless, a program of measurement of the power 
densities is likely to be a very long one. It may be desirable instead to estimate the power- 
density patterns from one of the theories of the scattering mechanism and a knowledge of the 
transmitting-antenna pattern. 
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Throughout the remainder of this section we will suppose that the power densities are to 
be measured. Questions as to when and for how long to average are left to the judgment of the 
proprietor of tlie receiving site; our object will be to suggest methods of measurement. 

One of the rectangular power densities, which we will call 




P. U,m) P,* (i + X, 


(33) 


n be expressed in terms of the arriving fields as follows: 


(f, m) = P. U, m) E * (0, 0, 0) 




c,y) dx dy 


<(>aik(^-y) = Ei^(x,y,0) E«^(0,0,0) 


(35) 


^ [We have used (Zl) and (2Z).] A method of measurement suggests itself immediately; Voltages 
proportional to the fields at two points on the xy-plane would be available at the terminals of a 
pair of appropriately placed dipoles. The Fourier transform of (35), the temporal cross- 
correlation of the fields, is the rectangular power density. The cross-correlation should be 
measured with one dipole located at the origin and the other located, in turn, at each of a suf- 
ficient number of points (x,y, 0) on the xy-plane to determine the behavior of through- 

out the plane. 



Figure 11 is a block diagram showing the essential parts of the circuitry that might be used 
to cross-correlate the dipole voltages V^ and V^; the simpler circuit of Fig. 12 could be used to 
obtain <#>^^(0, 0). the "cross-correlation" for identicaUy oriented dipoles at the same point in 
space. A radian carrier frequency <0 q is assumed for the arriving signal, and a local oscillator 
of radian frequency < co^ is shown in Fig. 11. It is assumed that the output circuits of the 
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PEAK DETECTOR 




Fig. 12. The a 


relator is just a watt-hour meter. 


four multiplying elements shown in Fig. 11 include filters that reject all but the difference - 
frequency component of the inputs. The balanced-modulator outputs are the direct ("d-c") 

ReiV^V^*} and tSmjV^V^*} , 

which could be averaged, for example, by Miller integrators^^ and recorded at appropriate 
intervals. In terms of these averages, the cross-correlation of the fields is 

</) ^.^(x, y) cc Re{ + j s®m{ 3^} , 

the proportionality constant being determined by the gain of the system. Again in Fig. 12 it is 
assumed that the output circuit of the squarer includes a filter to reject all but the second har- 
monic. The peak detector output !V^|^ should be similarly averaged to obtain a quantity pro- 
portional to 0). 

Only four rectangular power densities need be measured. For purposes of illustration, 
we will work with the four rectangular densities that would be measured with x- and y-oriented 
dipoles, although this choice is quite arbitrary. The "horizontal- and vertical-power densities," 
(32), can then be obtained from the measured relations (10). The nec- 

essary algebra, however, is somewhat tedious. At the risk of clouding the main issues, we will 
introduce matrix notation to simplify the algebraic manipulations. Suppose then, that we have 
the elements of the matrix 


^axx' ’ '^axy' ’ 

tl> u, m) _ ip u, m) 


of rectangular power densities, measured with x- and y-oriented dipoles only. Define the fol- 
lowing column matrices: 
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transformation matrix that has the coefficients 

The relation between P] and F] is given by th 
of the first two of equations (10) as elements; 



Note that the determinant of the T-matrix is 
I [T(f,m)]| = |- , 

so that [T] may be inverted freely except possibly for the directions f ^ + m^ = 1, corresponding 
to the horizon. According to (10), the relation between P] and F] is 

P1 = [T]XF1 ; F] = (T]-ixP] . <39) 

We will write the transpose of F] as 
F = P X [T] , 

reserving the star to mean that the complex conjugate is to be talcen. 

The matrix of power densities, (32), may now be written 

[*(i,m)) = ^ yy F^(i,m)]X F*(i + X, m + pi dX dp , (31a) 

"n [T(i,m)]"^ X P^(i,m)] X P*{1 + X. m + p) X [T(i + X, m + p)]^'^ dX dp . 

The elements of [T] will remain essentially constant over the range of appreciable contribution 
to the integral, so the T-matrices may be removed from the integral. We recognize the re- 
maining integral as (1, m)], giving the result 

[4(f,m)] = ^ [T(i,m)l"^X [;),^(f,m)]X [T(f,m)l^'^ . (40) 

To recapitulate, the cross-correlations of x- and y-fields on the xy-plane (35) are measured, 
these four quantities are transformed (34), giving the elements of the ip -matrix (36), which, 
when operated on by the T-matrix (38), gives the matrix of power densities (40). 

But there is a fundamental difficulty with this procedure for measuring the power densities: 
We have tacitly assumed that the total field at each of the measuring dipoles is the incident field 
arriving from the scattering region; no fields reflected from below the xy-plane are assumed 
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present. In other v.'ords, we have assumed that the region below the xy-plane is free space. In 
some cases, this assumption may be reasonable, but usually an antenna will be erected over a 
more or less level ground, which must be taken into consideration in measuring the arriving 
power densities. 


If the ground at the receiving site is smooth enough, if irregularities in the terrain are 
much smaller than a wavelength, the reflected fields are easily calculated. In this case, the 
foregoing measurement procedure may be re-interpreted to account for the presence of the 

ground. We will suppose that the xy-plane is 
parallel to the ground and at a distance h wave- 
lengths above the ground. Cross-correlations 
may be measured as before on the xy-plane or 
on some other plane that is tilted with respect 
to the ground. We will start by interpreting 
cross-correlation measurements on the xy- 
plane over a lossy, but otherwise smooth, 
homogeneous ground. 

Receiving-site geometry is shown in 
Fig. 13. It may be seen that the reflected 
waves arriving at the point (x, y, z) above the 
ground travel a greater distance than the cor- 
responding incident waves. This extra dis- 
tance is just accounted for by calculating the 
phase of a reflected wave as though it were 
not reflected, but continued instead, in the 
absence of the ground, to the image of the 
point (x, y, z). 

On being reflected, the wave will suffer 
an additional phase shift and attenuation, which 
will depend on the angle of incidence and the 
polarization. This effect may be accounted 
for by assigning reflection coefficients 
Pjj(i, m) and py(f, m) to the reflected horizontal and vertical plane-wave components. Sign con- 
ventions for the reflection coefficients are shown in Fig. 14. Note that the reflected horizontal i 

component PjjFjj^ contributes to the x- and y-fields with the same sign as the incident horizontal i 

component; the reflected vertical component py Fy^, however, contributes to the tangential fields h 

with the opposite sign from the incident vertical component. If the ground were a perfect con- i 

ductor, py = — 1 and py = + 1, so that both reflected components would contribute to the tangen- i • 

tial fields with opposite signs from the incident components. ! '• 

w 
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r [T(i, m)]“^ X [I/* U, m)l X [TU, m)]^ 


Given cross-correlation measurements on the xy-plane but no information about the sur- 
roundings of the plane or its orientation with respect to the scattering region, we are able to 
determine the direction cosines 1 and m of the power density impinging on the plane, but we 
are not able to determine the sign of the direction cosine n of the arriving power density. In 
other words, we are not able to distinguish between power density arriving from a point above 
the measurement plane and power density arriving from the image of that point below the meas- 
urement plane. We were able to interpret measurements on a plane parallel to the ground be- 
cause we knew that the scattering region must be above the plane and that plane waves arriving 
from below the plane must be reflected copies of the incident plane waves. If the measurement 
plane were tilted with respect to the ground, however, it would be possible for incident radia- 
tion to arrive from both sides of the plane. Nevertheless, it will generally be possible to orient 
a tilted measurement plane with respect to the scattering region in such a way that we are rea- 
sonably sure that the incident plane waves arrive only from one side of the plane. 

Next we will interpret cross-correlation measurements on a vertical plane above the ground; 
the xz-plane of Fig. 13. We suppose that the positive y-axis is oriented in the general direction 
of the transmitter, so that all incident radiation arrives with a positive direction cosine m. We 
will obtain a two-lobed power-density pattern: One lobe will extend above the ground in front of 
the xz-plane, and the other lobe, below the ground, will be the somewhat attenuated image of the 
first because of ground reflection. At first glance, it would appear that these two lobes, taken 
together, could be delivered to the antenna designer as the power-density pattern for the design 
of an optimum antenna to be erected at the origin of coordinates in Fig. 13. This, however, 
would be a mistake! A plane wave arriving from above the xy-plane is perfectly correlated 
with the corresponding plane wave that is reflected from the ground and arrives from the image 
direction below the xy-plane, so that the assumption of incoherence, discussed in Sec. II, is 
not valid for the superposition of incident and reflected plane waves. As we shall see, the two 
lobes will be scalloped by alternate constructive and destructive interference between incident 
and reflected waves. The angular locations of the scallops will depend critically on the height h 
of the reference dipole used in the measurement. In contrast, if any two arriving plane -wave 
components were uncorrelated unless they arrived from nearly coincident directions, as re- 
quired by the assumption of incoherence, the power densities would be, for all practical pur- 
poses, independent of the location of the reference dipole used in measuring them. 

The proprietor of the receiving site should measure the densities of incident power and de- 
liver these, together with the ground-reflection coefficients, to the antenna designer for his 
consideration. 

Returning to the problem at hand, we may write the column matrix of x- and y-field com- 
ponents on the xz-plane, which is normal to the ground as 
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= jy [TU,m)] x{[l]+e-j^’^<^+^^"[P(f,m)]}xF^U,m)] df dm (48) 

from (38), (41) and (44). The cross-correlation matrix to be measured is 


[<(>^{x,z)] = E^(x,z)]XE*(0,0) . <49) 

In writing an expression for this matrix in terms of the power densities, it will be convenient 
to convert one of the ingredients (48) into an integral over the direction cosines t and n, m being 
determined as the positive square root of 1 -i -n : 

di dm = ndS2 = ^ di dn 
In this v/ay, (49) becomes 

[,^(x,z)]= yy [T(i,m)]x{[l]4e-^^'(^^")"[p(i,m)]} ^ d. dn 

X yy F^(i,m)]XF*(i',m“) X{[1] + [p^(i'.m-)]} X [T(f',m')]^ di- dm' . 

As before, if the height h of the reference dipole is not too great, the last two factors may be 
removed from the second integral, which is then recognized as n[$]. The result, 

X [T] X {[1] + [p]} X [$] X{[1] + [p*l}X [T]^ . 

however, is not in the desirable form of a Fourier transform because of the function of z, 
g-j4ir(z+h)n^ appearing in the integrand. This difficulty may be remedied by splitting into 
the sum of two integrals: 

[<.n(x, Z)1 = yy ^ didnx [T] X [#] X { [1] + [p-*-]} X [T]^ 

+yy ^ ej2’r(xl-zn) ^ ^-j4xhn x [p] X [#] x { [1] + [p*]} X [T]t 
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itten in terms of its Fourier transform: 


r,/,n)].[.-(n)]} 


[:;,+(n)] = ^ [T(n)lX [$(n)]X{[l] H 


^ [T(-n)] X [p(-n)] X [$(-n)] 

X { [1] + [p*(-n)]} X [T(-n)]^ j 

se two , -matrices describe the two lobes mentioned in the preliminary discussion. The 
sence of the matrix, 

{[1] + [p*(n)]} . 

a,, sclloplng of to .PP.O lob. .no to coro.spoodln, lb tb. .bpr=..io~ lor Ip'l 

i<!PS the scalloping of the lower lobe. . 

It is important to note that, because of the presence of the ground, the matmx of arrmng 
.er densities. [*(n)]. will be zero for directions n < 0. below the horizon. Similarly. 

_nH Will be zero for n > 0. In other words, the upper lobe is zero below the >^orrzon. and 
: lower lobe is zero above the horizon. These facts allow the two lobes to be separated. T , 
,m the measured matrix of cross-correlations of x- and y-fields on the xz-plane. (49). we 
ly obtain the upper lobe: 




According to (50). the desired matrix of arriving power densities may then be calculated from 
j*] = J|. [Tf^ X [,(i'^]X [T]^‘^ X{[1] + . (5: 


UNCLASSIFIED 


\ 


/ 


UNCLASSIFIED 


If the distance between transmitting and receiving sites is very great, the signals received 
on the measurement dipoles may be so weak and corrupted by noise that it may be desirable to 
make the measurements instead with more directive antennas to improve the signal -to -noise 
ratio. In this connection, it should be noted that the long averaging times contemplated for the 
cross-correlation of antenna voltages will not be effective in removing that part of the "noise" 
in the measured cross-correlation which is due to coherent noises received on the two meas- 
uring antennas. Suppose that the input to the circuit of Fig. 11 from the left-hand antenna is 
+ (Nj + k^N^) 

and the input from the right-hand antenna is 
^2 + ■ 

and are the two signal voltages. (N^ + k^N^) and (N^ + k^N^) are the two noise voltages, 
and k^ and k^ are constants. With the circuit of Fig. 11 and associated integrators we will 
measure 


{V^ + (N^ + k^N^)} {V* + (N* + k*N*)} = V« + (N^ + k^N^) (N* + k* N*) , 

since the signals and noises will be uncorrelated. The noise components. N^. and are 
assumed uncorrelated; is that part of the noise which is common to both antennas. Note 
that the measured cross-correlation, which reduces to 

VTV* + k. k* In P , 


e contemplate using a circuit like that 


includes a term due to this common noise. 

Another way of looking at the noise problem is that v 
of Fig. 11 to measure the density of power arriving at the receiving site. In addition to the sig- 
nal power, we expect a more or less uniform distribution of noise power density, as discussed 
at the end of Sec. IV. Since signal and noise are incoherent, their power densities add, and 
we should expect to measure the sum of these two power densities with equipment appropriate 
for the measurement of power density. The advantage in using directive antennas for the meas- 
urement is that they would capture less of the total noise power, but essentially all of the signal 
power. To a first and fairly good approximation, a constant power density may be subtracted 
from the measured self-densities, 

The measurement may be improved considerably by transmitting an unmodulated carrier 
and using sharply tuned circuits to remove all but a narrow band of noise. The tuned circuits 
should precede the multiplying elements in Figs, ii and 12. 

It may also be desirable to make measurements with directive antennas for another reason: 
The ground at the receiving site may be so irregular that a calculation of the fields reflected 
from it would be impractical. This difficulty may be overcome by measuring the power densi- 
ties with antennas that are directive enough to reject, substantially, the radiation reflected 
from the ground. As before, a measurement of cross-correlations of antenna terminal voltages 
will be required; one antenna will be located at the origin of coordinates and another should be 
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located in turn at a number of points on the xy-plane. Since we will assume that neither antenna 
is sensitive to fields reflected from the ground, considerable freedom in the choice of orienta- 
tion of the measurement plane (the xy-plane) is left to the proprietor of the receiving site. This 
freedom will be restricted only by the requirement that all the incident radiation should arrive 
from above the xy-plane: with positive direction cosine n. 

In Sec. IV, we developed an expression, (29), for the open-circuit voltage at the terminals 
of a receiving antenna due to an arriving scatter signal. The reference point on the antenna to 


of a receiving antenna due to an arrr 
which the phase of its pattern is refe 
If we had assumed, instead, that the 
point on the antenna is at (x, y, z), we 


assumed to be located at the origin of coordinates, 
s located at the point (x, y, z), i.e., the reference 


) = ^ 2 y F.(i, 


The phases of the arriving plane-wave components are shifted by the same factor as in (21), 
the expression for the arriving fields. Note that, in moving this antenna from the origin to its 
new location, we keep it pointed in the same direction, and we do not rotate it, so that its polar- 
ization remains the same. 

In the dipole measurements we used x- and y-oriented dipoles; in the present measurements 
with directive antennas we will also need two kinds of antennas, which will have terminal voltages 
Vj(x, y, z) and V2(x, y, z), given by (53). Antenna 1 might be a small dish with horizontally polar- 
ized feed, and antenna 2 might be the same dish with a vertically polarized feed, for example. 

In place of the T-matrix, we will need the square matrix, 

Fjj^(f, m) Fh2^''’ 

[F(f,m)] = , (54) 

F„.(i, m) F„,(f,m) 


of horizontal and vertical plane-wave components emitted by the two antennas when they are 
excited with unit current. When the antennas are located on the xy-plane, the column matrix 
of their open-circuit voltages is 

V^(x,y, 0)1 

V(x.y)l = y [F(i,m)l^X F^(f,m)] dn , (55) 

V2(x,y,0) u.h. 
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[<p{x.y)] = V(x,y)]X V*(0,0) 


y [F(f,m)]^ej2’^(=^+y”' dfi 


F (f, m)] X F *(f', m') X [F*(f', m’)] dil' 


The range of both integrals n 
removed from the second inti 


; extended, and the matrix of antenna patterns, [F*], may be 
which is then recognized as [$(f, m)]. The result of these 


(>(x,y)] = fl [«f,m)]ej^ 


m. m)] = (^') ^ [F(je, m)] X [$(i, m)] X [F*(f, m)] 


The matrix of measured cross-correlations, (56), is transformed, (57), and operated on by the 
matrix of antenna patterns, (54), to give the matrix of power densities: 


[$(f, m)] = n(-^) [F(f, m)]^'^ X m. m)] X [F*(f, m)]' 


We have suggested a number of methods for measuring the power densities, all of which 
involve measuring the cross-correlations of antenna voltages on some plane, transforming, 
and performing algebraic operations on the transforms. It is time to take stock of the results 
to see which parts of these procedures will present real difficulties and, if possible, to suggest 
short cuts. The algebraic operations, being relatively simple and straightforward, may be 
dismissed immediately. 

The measurement of the cross-correlations will be the most difficult operation. Making a 
measurement at each point on the measurement plane, leaving the movable antenna fixed at each 
point for enough time to take an average, would be physically impossible. Even measuring at 
a sufficient number of points to determine the behavior of the cross-correlation throughout the 
plane will be a time-consuming operation. It will be advantageous to reduce the number of 
measurement points to the bare minimum. 
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With the exception of the measurement on a plane perpendicular to the ground, which will 
be discussed separately, all the results, (40), (47) and (58), have been of the form 

[$(i, m)] = [A(i,m)] X [^'(i,m)I X [A*(f,m)]j , 

in which [$] is the matrix of arriving power densities, [^1 is the matrix of Uie transforms of 
the cross-correlations — a matrix of rectangular power densities - and [A] is the appropriate 
transformation matrix. [#] is the matrix of a non-negative definite Hermitian form, as we 
discovered at the beginning of this section. Because of the symmetry of its expression in 
terms of [:p], it follows that [i/] is also the matrix of a non-negative definite Hermitian form. 
In particular. 


e the Fourier transforms of the 


tiij,(i,m) = ^*.(i,m) , 

so that the are real. The cross-correlations, which £ 

Ip's, must therefore obey 

<,ik(x,y) =0^.(-x,-y) , 

SO that the self terms are conjugate symmetric: 

<f'ii(x,y) = <;)*(-x,-y) . 

By making use of these symmetries, it will be possible to obtain the four power densities from 
a measurement of the four cross-correlations over only half of tlie xy-plane — an appreciable 
saving in effort. 

This saving is not possible for the measurement on a plane normal to the ground as we 
have developed it, because the matrix of transforms of the cross-correlations, [i)'''(n)] + [v"(n)], 
is not Hermitian, as can be seen from (50). The extra penalty attached to this measurement 
may be attributed to the fact that it is also a measurement of the ground-reflection coefficients: 
If [*] is eliminated from the equations (50), there results 

[p(n)] = [T(n)]'^ X [^■(-n)] X X [T(n)] 

Although this feature may be desirable, one could certainly devise simpler experiments for the 
measurement of the reflection coefficients. 

In some cases, it may be reasonable to assume another kind of symmetry among the ip's. 
due to symmetry about the plane containing the great circle path from transmitter to receiving 
site. If the y-axis lies in this great-circle plane, the cross-correlations will show some kind 
of symmetry about the y-axis, depending on the polarizations of both the transmitting antenna 
and the measuring antennas. In these cases, which wiU not be enumerated, it wUl be possible 
to obtain the power densities from a measurement of the cross-correlations over a quarter of 
^e measurement plane except where the measurement is on a plane perpendicular to the ground. 
In that case, the yz-plane should coincide with the great-circle plane, the symmetries will be 
about the z-axis, and a measurement over half of the xz-plane will suffice 

It is a f^damental property of Fourier transforms that the fine structure of a function is 
determmed by the gross behavior - the extent - of its transform, and conversely. Thus, the 
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required extent on the measurement plane of the cross-correlation measurements would depend 
on the amount of fine structure in the power densities if the object of the experiment were to 
determine the power de nsities as accurately as possible. But if the measured power densities 
are to be used in the design of an antenna, we should like to know how much of the fine structure 
is necessary for the design; perhaps some could be neglected with a resulting saving in meas- 
urement effort. As will be shown in the next sections, 
this question is best answered directly in terms of the 
necessary extent of the cross-correlation measurements 
on the measuring plane: It will be necessary to know the 
cross-correlation between arriving fields at any two 
points on and possibly slightly beyond the aperture of 
the antenna; if the antenna is to have a large aperture, 
the measurements should extend over a correspondingly 
large area. 

On the other hand, the fine structure of the cross- 
correlations, which will dictate the allowable separation 
of the measurement points, is determined by the angular 
extent of the power densities. If it may be assumed that 
all the plane -wave components of the signal arrive from 
within the rectangular range of directions AO, shown in 
Fig. 15, the familiar Sampling Theorem will suggest an 
appropriate spacing of measurement points and perform the Fourier transform of the cross- 
correlations measured at these points. 

One of the rectangular power densities, which will be zero outside of Afl, may be repre- 
sented by a Fourier series inside All. In fact, the product of this rectangular density with an 
arbitrary function of (f, m) may be represented by a Fourier series inside All: 

j2jr[x (f-i )+y (m-m )] -j27r[r'(i-f )/Af+s'(m-m )/Am] 

ip^^^U.m) e = Z Z ® 



Fig. 15. Measurement effort is reduced if 
the power densities arrive from within a 
narrow range of directions, AS2. 


so that the rectangular power density itself is 

-j2jr[(x +r'/Af)(f-i )+(y +s'/Am)(m-m )] 

<P^^U.m)= Y C^.g. e 

r', s' 

The coefficients C^,g, and the constants x^ and y^ are as yet undetermined. The cross-correlation 
</)^(x, y) is the transform of m): 

<l'ik<*-y) = yj dt dm 

All 

j2tr(xf^+ym^) sinir[(x - x^) Ai - r'] sinir[(y - y^) Am - s'] 

= AfAme 2 ^r's' ir[(x - x^) Af - r'] tr[(y - y^) Am - s'] 

r'. s' 
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d y = y + s/Am, all the terms in the double s 
: s. Accordingly, 


<(),, (x + r/Ai, y + s/Am) = 


j25r[(x +r/Ai)i +(y +s/Am)m 


j2;r[(x-x )i +(y-y )m ] 
(x,y) = e ° ° ° ° 2 


-j2jr(riyAi +smyAm) sm7r[(x - x^) Ai - 
' ir[(x - X ) Ai - r 


The assumption that the power de 
AS2 has made possible this expre 
terms of samples of the cross-cc 


e from within the rectangular 


urge of directions 
m the xy-plane in 


taken at a rectangular array of discrete points on the xy-plane. Rows of points are separated 
by l/Ai in the x-direction and l/Am in the y-direction. The smaller the range of directions 
ASJ of the arriving power densities, the greater is the allowable spacing of the measurement 
points. Even if we could only state with certainty that none of the storage fields of the source 
reach the receiving site, so that the power densities might arrive from any direction within 
+ m^ = 1, we could still inscribe that range in a square of sides Ai = Am = 2, and we could 
determine the cross-correlations everywhere on the xy-plane from measurements taken at a 
rectangular array of points spaced by a half wavelength. 

The constants x and y are still quite arbitrary (they should be real numbers); they may 
be chosen at the discretion of the experimenter to displace the measurement grid with respect 
to the origin, where the reference antenna is located, so that none of the measurement points 
coincides with the origin. 

The rectcingular power density, which is the Fourier transform of the cross-correlation, 
is just the Fourier series postulated at the start, which, in terms of the coefficients we have 
determined, becomes 

f 1 _ -j2jr[(x +r/Ai)i+(y +s/Am)m] 


for (f , m) inside A iJ 


Unfortunately, this application of the Sampling Theorem only applies directly to the meas- 
sment on a plane perpendicular to the ground if the range of directions Afl is large enough to 
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include the directions of arrival of both incident and reflected plane waves. This would tend to 
reduce the permissible vertical separation of measurement points on the xz-plane. 


The methods of power-density measurement discussed so far have had the common advan- 
tage that relatively modest antennas may be used. As one of these is shifted around on the 
measurement plane, the field pattern that would appear over the aperture of a larger antenna 
is explored. But the most intuitive device for measuring power density — a very directive (and 
hence very large) steerable antenna - has not yet been discussed. 

It should be clear that if a single, large antenna is steered across the sky and the received 
power measured at appropriate directions, only a single power density can result. With two 
such antennas, used separately, it might be possible to measure the self-densities and 

^YY‘ other two, cross-correlation measurements will again be needed. One attack 

is suggested by the first part of (34), which reads 

A very directive "i-polarized" antenna pointed in the direction (f, m) would have a terminal 
voltage proportional to P^^(jf, m), and a k-oriented dipole located at the origin would have a 
terminal voltage proportional to Ej^^(0, 0, 0). Cross-correlating these two, we would obtain 
one of the rectangular power densities. Unfortunately, there are a number of difficulties with 
an experiment of this type. One of these is that, when a practical steerable antenna of finite 
aperture dimensions is used, it is advantageous to mount the steerable antenna close to the 
reference dipole - so close, in fact, that the structure of the steerable antenna is just in front 
of or just behind the reference dipole, interfering with its reception. 

A more convenient measurement can be made with a single, large, steerable dish having 
both horizontal and vertical feeds at the focal point. The feeds should be arranged so that the 
mutual impedance between the horizontal terminals and the vertical terminals is zero. This 
allows the horizontal terminals and vertical terminals to be regarded as the terminals of two 
separate antennas, each steered in the same direction and located at the same point in space 
(the origin of coordinates). These antennas may be described by the square matrix, 

■FHh(^.m,^,e) F^^(i,m,^,e)' 

[Fg(i,m,<p,e)] = 

F.y.jj(i,m,^,9) 

of horizontal and vertical plane-wave components emitted when the antennas are excited with 
unit current. Fjjj^ and F^^^ are the horizontal and vertical patterns resulting from excitation 
of the horizontal feed; F^^ and Fy^ are the horizontal and vertical patterns resulting from 
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excitation of the verticai feed. To a good approximation, and B are zero, but we sha 

retain them for completeness. As has been indicated, the elements of the pattern matrix will 
depend on the angles <p and 0 through which the antenna structure has been elevated and rotated 
in azimuth. 

The open-circuit voltages of the two sections of the antenna, according to (55), are 


V^(^P,e)1 

V^(<p,e) 



u.h. 


FJf,m)]dn . 


This time tlie correlation matrix to be measured is 

e)] = Vg(<p,e)ix V|((p,e) 

= (^f y [F^(i,m.<,,e)]^dn 

X ^ F^(f,m)] X F*(i',m') X [F|(i', m’ , (p, 0)] dfl' , 

uJi. ' ' 

and, as before, the pattern matrix may be removed from the second integral.which is then rec- 
ognized as [*(f.m)l. The remaining expression, 

0)1 = (^)^ y [Fg(f,m, <p,e)l^X[*(i,m)]X [F|(f,m, <p,e)] dJ2 . 

will have to be treated somewhat differently this time, however. It will be assumed that, when 
the antenna structure is elevated by an angle <p and rotated by an angle 0, the patterns have sharp 
maxima in the direction ~ power-density patterns remain essentially 

constant over the range of directions for which there is an appreciabie contribution to the 
integral. The integral becomes 

Us(<P.e)] = y [Fg(f,m,,p,e)]^X[*(f^,m^)]X[F*(i,m,».,0)]dR . (62) 

u. h. 

In this form, the ♦-matrix is constant so far as the integration is concerned. Unfortunately, 
it cannot be removed from the integrand because of the variable matrices that pre- and post- 
multiply it. Nevertheless, the 16 cross-products that multiply the various elements of the 
♦ -matrix may be integrated, and from these and the measured elements of the correlation 
matrix [#g] the four power densiUes may be obtained. If F.^jj and Fjj^ may be neglected, only 
four integrals, one of which is the conjugate of another, need be evaluated; each power density 
is then proportional, through one of the four integrals, to the corresponding correlation. 
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From (62), the density of power arriving from (i^, m^) may be obtained. The measurement 
and calculations should be repeated for other directions until there are enough data to determine 
the power-density patterns in their entirety. Strictly speaking, the pattern integrals in (62) 
should be re-evaluated each time the elevation angle tp is changed. However, if the patterns 
are extremely peaked and if no measurements are to be made in the general direction of the 
zenith, this may not be necessary. 

The same kind of long-time average measurement must be made with the antenna pointed 
in one direction as is necessary when the movable antenna is located at one point in the alter- 
nate technique of making cross-correlation measurements on a plane. And, in general, the 
number of d irections in which measurements are made should agree with the number of meas- 
urement points chosen for the alternate scheme. Thus, although the use of the very directive, 
steerable antenna for power-density measurement is conceptually more straightforward than 
the cross-correlation measurements on a plane, more elaborate equipment is required, and 
the effort involved will be about the same. 
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VI. OPTIMIZATION OF AN ARRAY 

The design of two specific types of receiving antenna will be considered — one in this section 
and another in Sec. VII. In each case, the required information about the arriving scatter signal, 
which is given by the power-density patterns, is presumed known. The two antennas will differ 
because different restrictions will be placed on their construction. It is apparent that some re- 
striction should be placed on the shape, and particularly on the size, of the aperture: In general, 
the larger the aperture dimensions, the more power may be extracted from the arriving scatter 
fields, although, as we saw in Sec. I, the point of diminishing returns is reached when the aper- 
ture of the receiving antenna is as big as the transmitting-antenna aperture. 

The antenna to be considered in this section is an array of smaller, not necessarily identical, 
antennas. The number, construction, location and orientation of the array elements is presumed 
to be included in the list of restrictions on the construction of the antenna; although we will not 
specify these choices, they are important and 
should be given due consideration. The re- 
mainder of the design, which will be considered 
here, will amount to the design of the inter- 
connections between the terminals of the array 
elements and the common load to which the re- 
ceived power is to be supplied — a network 
problem. 

As shown in Fig. 16, we will start by de- 
signing a linear coupling network that will ex- 
tract the maximum mean power from a group 
of disconnected sources and supply it to a single Fig. 16. The coupling network for a group of discon- 
load. N sources, each having an open-circuit nected sources. 

voltage V., an internal resistance Rj, and supplying a current I., are indicated on the left of the 
coupling network. The single load to which the power is to be supplied is R^ on the right; I^ is 
the load current. The extra voltage source V^ on the right does not belong in the description of 
a receiving antenna; it will be used temporarily to help demonstrate a restriction on the elements 
of the coupling network. As usual, the V. and the L are the time-varying complex amplitudes of 
the corresponding voltages and currents. In practice, the coupling network might consist of a 
collection of transmission lines that pass all the frequency components of a typical communica- 
tion signal without appreciable relative attenuation or phase shift: We will assume that the input- 
output relations of the coupling network, calculated for steady-state sinusoidal inputs, will hold 
when the inputs arc fading, modulated carriers. 

The circuit of Fig. 16 is a fairly accurate representation of the interconnections between the 
elements of an array and their common load. V. may be regarded as the open-circuit voltage of 
the i*’’ element and R. its radiation resistance, so that the i**^ source on the left is the Thevenin 
equivalent of the i*^ element as seen from its terminals. (It is assumed that the driving-point 
reactance of each element has been resonated.) But if two antennas are located close together 
as they might be in an array, and if one is driven with a generator, a fraction of the generator 
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tj = 2 Y^. (65) 

is the transmission coefficient from the i^^ source on the left in Fig. 16 to the load (or from the 
right-hand source to the i**' source resistance on the left). We have just demonstrated that 
N 

E , (66) 


which is equivalent to the statement that the sum of the squared magnitudes of the elements in 
the first row of the scattering matrix^^ is less than or equal to one when (in the present case) 
the load reflection coefficient — the diagonal element — is zero. 

Evidently the inequality (66) becomes an equality if the coupling network is lossless. It will 
also be an equality if the coupling network contains resistors that are not coupled to the load , 
but that may, however, be coupled to one or more of the left-hand sources. 

We will have no further use for the right-hand voltage source, which may be considered to 
be removed (V =0). The mean power in the load from the left-hand sources is 
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would only be the fraction of that possible if (66) were satisfied with the equal sign. Thus, 
we should adjust the transmission coefficients in such a way as to maximize (67) subject to the 
constraint, 

N 

S |t. |^= 1 . (68) 

We have shown that the constraint (68) is a necessary condition on the realizability of the 
optimum coupling network, but as yet we have not shown that it is sufficient . To show that the 
constraint is sufficient, it will be simplest to produce a realization of the coupling network, 
which, of course, is one of the desired results. Such a realization is shown in Fig. 17. The 



Fig. 17. A realization of the optimum coupling network. 


boxes labeled are phase shifters that are coupled by ideal transformers to the load R^. 
Physically, the i phase shifter might be a transmission line of characteristic impedance ° 

R. and electrical length 0. radians. The i*** ideal transformer, which has a turns ratio of a. to 
1, might also be approximated by appropriate connections to a resonant length of transmission 
line. The parallel rather than the series connection of transformer secondaries has been chosen 
to simplify the construction of the coupling network: A single autotransformer (a single resonant 
line) with various primary taps will serve the purposes of all N ideal transformers shown in 
Fig. 17. 

With the load short-circuited, the current flowing in the short is 



v 
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Comparing (69) with (67) and (70) with (68), we recognize that the transmission coefficients for 
the network of Fig. 17 may be taken to be 



If the turns ratios and phase shifts are appropriately chosen, the network of Fig. 17 will indeed 
be a realization of the optimum coupling network. 

It remains to maximize the mean load power by adjusting the transformers and phase shift- 
ers. We will work with the transmission coefficients, from which the turns ratios and phase 
shifts may be determined by (71). Manipulating (67), the mean load power becomes 
N N ^ ^ 

W^=l/8 I E titk"-7== 2 2 ^iV-^ik • 
i=l k=l -s/^i^k i=l k=l 

And we are not surprised to discover that the data required for the maximization consists of the 



* - V V* 
8JR^ 


(73) 


of quantities proportional to the temporal cross-correlations of open-circuit source voltages. 

We presume that these quantities, which may be determined from the power densities by the 
methods of Sec. V (or which may be measured directly! ), are known. In any case, we note from 
(73) that so that the auto-correlations are real and non-negative. 

We will use the method of Lagrange to maximize (72) subject to the constraint (68): We will 
maximize 

N N N 

U = 2 ‘i 2 “ 2 |t. I 

i=l k=l i=l 
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where a is the Lagrangian multiplier. There are 2 N relations; N state that U is stationary with 
the N real parts of the t's, and N state that U is stationary with the imaginary parts. Tempo- 
rarily, we let + j cTj. where the r and are real. Then 

N N N 

U= Z + E + ' 

i=l k=l i=l 

and the first N relations are 

N N 

= Z ‘'’nk'^ 2 (^i + j<^i)^in-2“\ = 0 

k=l i=i 

for n = 1, 2 N 

The second N relations read 

N N 

1 (^k-j<^k> “^nk" j S (’■i + j<^i)<?>in-2acr^ = 0 

" k=l i=l 

for n = 1, 2, . . . ,N 

The sum 

" ” i=l 

utilizes all 2 N relations in N complex relations: 

N 

^ = ° • n=1.2,...,N . (74) 

i=l 

At this point, it will be convenient to convert to matrix notation. Let the column matrix of 
transmission coefficients and the column matrix of normalized open-circuit source voltages be 
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Evidently, the constraint may be written 

J_Xt*] = l ; (68a) 

the mean load power is 

X [()> ] X t*] ; (V2a) 

and the result of the maximization reads 

X [())] = a (74a) 

Post-multiplying both sides of (74a) by t*], we have 
^X[<).IXt*] = a^Xt*]= a , 

so that the Lagrange multiplier a is the maximum mean load power we seek. 

The solution of (74a) for the maximum multiplier and corresponding t-matrix is the well- 
known characteristic value problem. Only enough will be said about it to indicate the nature 
of its solution and to show the physical significance of certain special cases. 

The result, (74a), may be written 

^X{[.#>]-a [1]} =0 , 

in which [1] is the unit matrix of order N. This equation has non-trivial solutions only if the 
determinant 

[{[</']- a [11}| =0 . 

that is. 
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which is the characteristic equation of the matrix [(()]. Evidently, this equation has N a-roots 
or characteristic values. In spite of the fact that the off-diagonal elements are, in general, 
complex, (72a) is a non-negative definite Hermitian form, so the characteristic values must be 
real and non-negative, and the coefficients of the characteristic equation must be real, and, if 
none of them is zero, they must alternate in sign. We are interested in obtaining the ma.ximum 
power; the greatest characteristic value is that maximum power. 

If the array elements have been properly placed, will be a distinct (not multiple) root 
of the characteristic equation. In that case, the rank of the system determinant 


fill 


will be equal to (N - i), and at least one of the cofactors C.^ of its elements will be different 
from zero. (74a) will then be satisfied by the characteristic vector (the set of transmission co 
efficients) 

and the column k may be chosen so that at least one of the C., is different from zero. In other 

Tors Th of one column of the matrix of cofac- 

tors. The constraint, (68), will also be satisfied if the constant of proportionality is chosen so 

t ^ .10 


I s iCi,r 


all of them. transmission coefficients to within an arbitrary phase shift 0 common to 

system relrmiLnrwnrbe“”ess d T ‘he 

be indeterminate. I„ this c e 4a, ' T 

‘al - t,, or by any linear combLal:o:It: "" “ ^ 
tl=atj + by . 

from the load. For this to be disconnecting the i‘h source entirely 

"""Is iTextrr''' P' ‘h^ 

for a three-element array baX^th suppose that a coupling network is desired 

ay having the matrix of normalized cross-correlations. 
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The reduced t-matrix will be taken from the first column of this matrix of cofactors, so that 
the entire t-matrix is 


It may also be seen by inspection that 


4 1 + j -i2 
l-i 5 1 + j 

j2 i-j 4 


= 6,-1, -1 - j. 0, 


so that 6 watts are still available after the elimination of the third source. In normalizing the 
characteristic vector (78), the phase shift 0 = 3ir/4 will lead to the simplest network. The 
t-matrix is then 


1 ^-i^/4 
nTI 


Jf 


and the corresponding network, for assumed equal source and load resistances, is shown in 
Fig. 18. 


v,tC^ 


45® PHASE SHIFT 
(i/8 wivelength line) 






3 





1-33-55551 


Fig. 18. In this example, one source may be entirely disconnected 
from the load. 
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As the number of sources (the order of the associated correlation matrix) increases, the 
labor involved in the method of solution just demonstrated soon becomes prohibitive. Fortunately, 
an iterative procedure, which is particularly well-suited to the problem of finding the greatest 
characteristic value and corresponding vector, is available. Automatic computing machines 
may be programmed to follow this procedure. 

The case of perfect correlation between source voltages and the case of complete incoher- 
ence of these voltages are two extremes worth mentioning. As discussed in Sec. II, if and Vj^ 
are "perfectly correlated," they will "fade together" with fluctuations in the scattering region. 

We may write the column matrix of normalized open-circuit source voltages (75) as 


= g(t) • v) = g(t) 


''n 

where the v. are not time functions. The matrix of normalized cross-correlations (76) becomes 
[<#>]= |g(t) ■ V] Xjr" - v) . 

In this case, the characteristic value problem is particularly simple. (74a) reads 
I j_x[<).l = _^xv]x^= , 

but this time, the product ^ X v] is a scalar constant, so that the result of the maximization is 
t oc V*. As before, we may choose the phase of the proportionality constant at our convenience, 
and its magnitude should be chosen to satisfy the constraint, (68a). Thus, if 

XXv*l = W^ , (’9) 

then may be chosen as 

t - v! 

The maximum available power, (72a), is 

V] X v;;' X t*] 

= ^ )L*Xv)x^^*Xv] = W^ . 

The maximum available power, (79), is seen to be the sum of the powers available from each 
source separately. As may be seen by expanding the characteristic equation, (77), the coefficient 
of (-a)*^‘^^ which is also the sum of the characteristic values, is just the sum of the diagonal 
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elements: the sum of the powers available from each source separately. Since one character- 
istic value is equal to the latter sum, and since the characteristic values are non-negative, all 
the other characteristic values must be zero. The case of perfect correlation just considered 
might arise if a small array were illuminated by a point source or by a relatively concentrated 
source having a power-density pattern that is sharply peaked in one direction. It might also 
arise if the object of the design were to produce a coupling network that is optimum only for a 
period of time short compared with the fading period of the arriving signal. The latter possibil- 
ity, which leads to a time-varying coupling network, will be considered in Sec. Vlll. 

The other extreme case is that of complete incoherence of the source voltages. In this case, 
the cross-correlations between the source voltages are zero, and the correlation matrix re- 
duces to the diagonal form: 



The characteristic values of this matrix are just its diagonal elements. If the greatest of these 
is the maximum power of <#>.. watts may evidently be obtained by matching the i*^ source 
directly to the load and disconnecting the other sources entirely. Almost complete incoherence 
would be expected for an array of widely separated elements that is illuminated by a diffuse 
source having a broad power-density pattern. 

For example, if the array is receiving only noise, and if that noise arrives uniformly from 
all directions in space, we would expect the array element voltages to be completely incoherent. 
We discovered in Sec. IV that the available noise power is independent of the antenna pattern; 
the noise power must be independent of the connections in the coupling network for the array. 

Since the characteristic values of the correlation matrix represent powers available with various 
connections, all the characteristic values must be equal. As may be seen by imagining that the 
sources are coupled to the load two at a time, this implies that the off-diagonal elements of 
the correlation matrix are zero, and that the diagonal elements are equal. 

On the other hand, suppose that two of the elements of the array are short dipoles. In the 
last section, it was shown that the cross-correlation of their open-circuit voltages is propor- 
tional to the Fourier transform of a rectangular power density. Since noise power arrives only 
from the real directions, 1, this rectangular power density cannot be a constant: It 

must be zero outside off^ + m^ = 1 and different from zero inside. Its transform, the cross- 
correlation, therefore, cannot be an impulse in space that is zero unless the locations of the two 
dipoles coincide. In other words, the open-circuit noise voltages of the two di[>oles will generally 
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be partially correlated! This apparent paradox has resulted from our neglecting the mutual 
impedances between the array elements. When these are properly accounted for, the noise 
power is truly independent of the connections in the coupling network between the array elements 
and the load. 



Instead of representing the array, as seen at the terminals of its elements, by the set of N 
disconnected sources shown at the left in Fig. 16, we are led to the more general representation 
of Fig. 19. The self and mutual impedances between array elements are the elements of the 



matrix [Z] of open-circuit driving-point and transfer impedances measured at the terminals of 
the array elements before these have been coupled to the load. The open-circuit voltages of the 
array elements are the elements of the column matrix 



The open-circuit impedances and voltages have been chosen for convenience; in some cases, tl 
open-circuit voltage of a dipole-array element can be calculated under the assumption that the 
other dipole elements, which are open-circuited, are removed entirely. 


] 
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Again we are to design a coupling network that will extract maximum power from the array 
of coupled sources and deliver it to a single load. We will attack this problem by reducing it to 
the one we have just solved: the design of an optimum coupling network for a set of disconnected 
sources. The first step is to resonate the self and mutual reactances. If the resonating network 
is connected in "series" with the array elements as shown in Fig. 20, and if the self and mutual 
reactances of the resonating network are of opposite sign from the corresponding array react- 
ances, the self and mutual impedances measured at the new terminal pairs are the self and mu- 
tual resistances of the array. We are left with the resistive network on the right in Fig. 20. 

It should be emphasized that we are dealing strictly with terminal pair behavior: There 
would be no connection, for example, between terminal 1 and terminal i in the resonating net- 
work; current can flow only to terminal i' from terminal i. It may be imagined that terminals 
i - i' are the ends of the primary winding of an ideal transformer, and that all other connections 
in the resonating network are made to the secondary winding. In any case, it should be clear 
that the open-circuit voltages of the resistive network we have obtained are the same as the 
corresponding open-circuit voltages of the array. 

Note that if we started with the resistive network on the right in Fig. 20 and added a "series" 
reactive network described by the matrix +j[X], we would reproduce the representation of the 
array (Fig. 19). Thus, if we can find a coupling network that extracts maximum power from the 
resistive network and delivers it to a single load, the resonating network and that coupling net- 
work are an arrangement for extracting maximum power from the array, for, if more power 
were available from the array alone, the coupling network could include the "series" reactive 


network +j [X] to obtain that greater power. 

The logic of the preceding argument hinges on the fact that the process of adding the reso- 
nating network is reversible. We shall have to use the same type of argument again, because 
we have not discovered how to extract maximum power from the resistive network. 
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We will call the matrix of output resistances [D] in anticipation of its being diagonal: 

E] = [D] X I] ; [D] = [b] X [R] X [b]^ (84) 

Evidently, a set of transformer connections described by any real matrix [b] can be realized 
if b . = 0, the corresponding transformer should be removed, leaving its primary circuit open 
and^its secondary short. Because of reciprocity, the resistance matrix [R] is symmetric, and 
it may indeed be diagonalized by the operation (84) with considerable freedom in the choice of 
the elements of [b]. 

Guillemin^^ discusses a method for obtaining the diagonalizing matrix, [b], as the product 
of M = N(N- l)/2 elementary matrices: 

[b] = [bj^] X X ... X [b^] X [bj . 

The elementary matrix [b^] is chosen so that [b^] X [R] = [R'J has a zero m the second row, first 
column: 


0 0 1 . 


^11 ^12 ^13 ■ 


^21 ^22 ^23 ■ 
^31 ^32 ^^33 ■ 


^11 ^12 ^13 ■ 

0 R- 2 R' 3 . 

R3^. R32 R33 . 


inwhichR^j =^2i~Rf7^1i- [b2] is chosen so that [b^l X [R'] = [R") has a zero in the third 

row, first column: 


'•31 "■32 ^33 • 


\i ^12 ^13 ■ 


in which = R3. - ^ R^.. Next, [b3] is chosen so that [bj] X [R"] = [R'"] has a zero in the 
third row, second column: 
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in which R'^’. = R'^. - ^ R'2i ■ The process continues into the fourth and lower rows until the 
remaining matrix! [b] xVl/is triangular. The diagonal elements of this matrix (R^^, R'22. 
R-''3. . . ) are the elements of the diagonal matrix [D], the source resistances on the right in 
Fig. 21: 


[D] = [b] X [R] X [b]^ 


When this diagonalization procedure is employed, the determinant of each elementary matrix, 
and hence the determinant of [b), is unity, [b] has the triangular form 


.0 no morn th.n N(N- O/Z .r«»(onn..rs nod N dlnec. conn.o .on. 

of .r„.Iorm.r conn.cUon. d..nrlb.d by [b|. 1W= ..vins of .on.ewl.M o.nr >■ J- 

.ibl, .r.n.torn..r. 1. . derided .dv„..*e .o .he n.e,hod of d.M,onM„....on nnder d..on...o„. 


UNCLASSIFIED 


UNCLASSIFIED 



Since it represents a physically reeilizable network, the matrix [D] of the resistive network 
as seen at the output terminals of the transformers must be non-negative definite. In addition, 
since it would be unreasonable to expect infinite power to be available from the array, each 
source resistance must be greater than zero, or if = 0, then V! = 0. The case D. = 0 and 
V! = 0 would not be expected in a practical situation; it could arise, for example, if some of the 
array elements were interconnected through a lossless network in such a way as to create an 
extra pair of terminals. VVe may expect, then, that [D] and [R] are positive definite, and that the 
diagonal elements D. are all greater than zero. It will be noted that the condition D. ^ 0 for 

i = 1, 2 k is necessary and sufficient to allow the diagonalization procedure to be carried 

through the (k + row. 

Since the determinant |[b) | = 1, [b]"^ must exist. (It has the same triangular form as [b].) 

(8 3) and (84) may then be rearranged to read 

V} = [b]"* X V] 

[R] = (bf* X [D] X [b]^'^ _ (84a) 

These equations show that if the disconnected sources on the right in Fig. 21 are followed by the 
set of ideal transformer connections described by [b]'^, the terminal-pair behavior of the re- 
sistive network and its open-circuit voltages are reproduced. In other words, the process 
brought about by the set of transformer connections is reversible. Thus, the set of tra.nsformer 
connections, followed by a coupling network that extracts maximum power from the disconnected 
sources and delivers it to a single load, is an optimum arrangement for extracting power from 
the resistive network, for if more power were available from the resistive network, the coup- 
ling network could include the set of transformer connections described by (b] * to obtain that 
greater power. 

We have shown that if the array is followed by a resonating network (Fig. 20), then by a set 
of transformer connections (Fig. 21) and, finally, by a coupling network (Fig. 17), the combined 
circuit will deliver the maximum power from the array to a single load. If the signal data con- 
sist of the elements of the matrix of cross-correlations of array-element voltages, V] X V*, 
the matrix of cross-correlations of the open-circuit voltages of the disconnected sources, shown 
on the right in Fig. 21 , is 

V] xy*. = [b] XV] X^^ X [b]^ 

from (83). And the matrix of normalized cross-correlations to be used in the characteristic 
value problem associated with the disconnected sources is 

1<P 1 = 1/8 X [b] X V) x;^ X [b]^ X , (85) 

where the diagonal matrix 
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.^^D, 


/D2 




The <(, -matrix has again been normalized in such a way that its characteristic values are extreme 
available powers. As usual, the form of the final coupling network is dictated by the solution 
vector corresponding to the greatest characteristic value. 

It is somewhat disheartening to contemplate the hardware required by the foregoing metho 
of extracting maximum power from a large array. Consider a 100-element array consisting of 
10 rows and 10 columns of antennas. The number is 10,000; about 5000 transformers and a 
similar number of elements in the resonating network would be required. Any such design 
would clearly be uneconomical. Of course, the mutual impedance between widely separate 
array elements would be small, so that many of the mutual impedances could probably be neg- 
lected altogether. This would effect a considerable saving in both the number of transforme 
and the number of elements in the resonating network. The mutual impedances -uld be par- 
ticularly small for an array of moderately directive elements that were arranged as far as pos 
sible not to have overlapping storage fields. Nevertheless, a number of circuit elements greatly 
in excess of the number of array elements would probably be required in any casm 
The advantage of the foregoing development has been that we were able pro 
to the maximum available power. Fortunately, the circuits we 

that can collect maximum power from an array; there are many sampler ones^ In fac , .ny 
circuit that, when connected to the array and the load, produces the same load current, wU 
serve just as well. The load current, from (63), is 


Io= 2 


Y . V' 


k=l 

to the disconnected sources and their coupling network. 
The admittances and voltages pertain to the 

In terms of the transmission coefficients, (65), 


n matrix notation. 


° 2Jr; 1^1 
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Finally, using (83), 


° ^4% ^ 


< [b] X VI 


The load current has been expressed in terms of (the normalized solution vector of the char- 
acteristic value problem), the source resistances of the separate sources, the matrix of turns 
ratios that diagonalizes the resistance matrix of the array, and the open-circuit array element 
voltages. Any circuit that delivers a load current that is the same weighted sum of open-circuit 
array element voltages will deliver the same maximum power as the circuits we have developed. 
This new circuit must provide the same weighting coefficients (transfer admittances) as appear 
in (87). 

The calculation of the transfer adm.ittances for the new circuit can also be simplified; it is 
not necessary to go through the intermediate step of calculating the matrix [b] that diagonalizes 
the resistance matrix, as we shall now see. The characteristic value problem for the array 
consists of determining the set of transmission coefficients^ that satisfies X [<#>] = a^with 
the greatest characteristic value With [<(>I given by (85), the equation reads 


^-^Jx(b]X{iv]XtL*}X[b],x[^^J=.X 

tiplyinj 

X [b]| X V)XV^} = c 


If the right-hand side is post-multiplied by 

[1] = [R]'* X [R] = [b]^ X [D]’‘ X [b] X[R] 
from (84a), we have 


\ t X — ^ X [bll x{4 V] XV*} = a-lix X[bllx [R] 


which may be written 

y X{| VIXV^} = ay X [R] , 


r X [b] ; 1, = y X [b]’^ X ■ 


The constraint on the transmission coefficients, x t*] = 1, becomes 

y X |b|'* X [-tj-* X f-tj-* X (b| ■* X y.) . 1 , 
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and, with the aid of (84a), it may be written 
^X[R]Xy*l = l , 

The lead current, (87), written in terms of the transfer admittances, i 


I = 




XVI 


In designing an optimum coupling netv/ork for an array of coupled sources, therefore, the 
characteristic value problem, (88), should be solved for the characteristic vector y], corre- 
sponding to the maximum characteristic value a^, which is the maximum available power. The 
characteristic vector should be normalized in such a way that the constraint, (89), is satisfied; 
as before, this operation will fix the magnitudes of the y's, but leave a common phase shift un- 
determined. Finally, a coupling network must be found that will deliver the load current, (90), 

for any set V] of open-circuit array element voltages; the transfer admittances are the set 

The characteristic value problem, (88), is more general than the one we were originally 
led to. Nevertheless, if (88) is post-multiplied by [R] : 

2X{| V] X V^ X[R1‘H = < 

it becomes formally identical with (74a) and may be solved by the same methods. The matrix 
{— V] X V* X [R]'^} , however, is generally not Hermitian, and it may be more convenient to 
solve (88) directly for its maximum characteristic value and corresponding vector. The fact 
that VpTv* is non-negative definite Hermitian and [R] is positive definite symmetric is sufficient 
to insure that the characteristic values are real, finite, and non-negative. This more general 
characteristic value problem is also well known and may be solved by iteration. 

One of the simplest forms the coupling network may take is shown in Fig. 22. The box 
labeled 6. is a transmission line of electrical length 6. radians and having any convenient 
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characteristic impedance; it may be simply an extension of the feeder from the i‘'’ array element. 
The ideal transformers may be realized approximately by appropriate connections to a single 
length of resonant transmission line. And if this resonant line is cut long (or short, as required), 
it will also perform the function of the resonating element. The design of this circuit is straight- 
forward when the characteristic value problem has been solved, making the coefficients in (90) 
available . 


Before closing this section, we will choose a simple array and show that the power available 
from it, when it is subjected to noise that arrives uniformly from all directions, is independent 
of the connections in the coupling network that connects the arr ay to it s load. For this to be the 
case, the cross-correlations between the open-circuit voltages Vl^ VJ^* of the disconnected sources 
on the right in Fig. 21 must be zero, and the same noise power must be available from each source. 
In other words, 

[■#>]“ [1] • 

Alternatively, as may be seen from (88), it must be that 
V] oc [R] , 

which we will demonstrate. 

As an array, we will choose a pair of short vertical half-dipoles just above a perfectly con- 
ducting horizontal plane and driven between their lower ends and the plane. The dipoles will be 
illuminated with noise which arrives uniformly from all directions in the upper hemisphere. If 
the dipoles have effective heights of a^ and a^ meters, and if they are separated by r wavelengths, 
the matrix of their open-circuit driving-point and transfer resistances is 


[R1 = ^ 


where 


f(r) = 


sin 2irr 
Zjrr 


cos 27rr _ sin 2irr 
(2irr)^ (2jrr)^ 


f(o) = lim f(r) = 2/3 
r-*0 


The open-circuit voltage of one of the dipoles is proportional to the z-component of the elec- 
tric field at the dipole. The cross-correlation of open-circuit voltages is therefore proportional 
to the cross-correlation of z-fields, which, as is suggested by (34), may be calculated as the 
Fourier transform of the rectangular power density From (33) and (10). 
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n 


'''azz*'*-'^' = ^ Fva‘''“> ^Va<^ + \, m + p) dX dp . 

The integral, as may be seen by comparing it with (31), is just $yy(f,m), the density of arriving 
vertically polarized noise power. As we discovered in Sec. IV, this quantity is a constant for the 
real directions in the case of uniform noise illumination; we shall again call it 



The cross-correlation of z-fields is 

- ff ^azz(''™> df dm 

= df dm 


The integral is again best evaluated in spherical coordinates with the aid of Watson. With 
I = p cose, m = p sine, it becomes 
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Allowing for reflections from the conducting plane, the matrix of cross-correlations of open- 
circuit voltages of the dipoles is 


V] X Y!l! = 


f(0) a^a^ f(r) 

a^a^ f(r) a^^ f(0) 


which is indeed proportional to their resistance matrix, as expected. We have shown that the 
noise power available from this particular array is independent of the connections in the (loss- 
less) coupling network connecting it to the load. 

Incidentally, the foregoing argument suggests an unusual way of calculating the real part of 
the mutual impedance between two antennas. 


78 

UNCLASSIFIED 


UNCLASSIFIED 


VII. OPTIMIZATION OF A CONTINUOUS ANTENNA 

In this section, we will discover how to adjust a receiving antenna of the type described in 
Sec. Ill in such a way that maximum signal power is available at its terminals. Instead of ad- 
justing a discrete set of transmission coefficients, as was done in Sec. VI, we will allow the 
adjustment of the electric field produced at each point on the aperture plane of the antenna when 
a unit current is applied to its terminals. In contrast to the array considered in Sec. VI, the 
antenna considered here is described by a continuous field distribution; it may be called a con- 
tinuous antenna. 

There is a pitfall in the path of the present investigation that we will dispense with at the 
start. Even if the region of space containing the active conductors of the antenna were speci- 
fied — in fact, even if this region had maximum dimensions small compared with a wavelength — 
it would be generally possible to fit i.nto this region an antenna having an arbitrarily sharp pat- 
tern! An antenna having a radian beamwidth small compared with the reciprocal of its largest 
aperture dimension, expressed in wavelengths, (or hav- 
ing a small beamwidth if the largest aperture dimension 
is small) is known as a superdirective antenna. Such 
an antenna is undesirable because its adjustment is crit- 
ical and its losses are relatively high. We will neglect 
antenna losses, and we will assume that the optimum 
configuration we arrive at analytically can be duplicated 
exactly in practice. Since, as we discovered in Sec. IV, 
the greatest signal power is available from the antenna 
having the greatest gain in the direction of maximum 
power density, we might expect our maximization of 
available signal power to lead to a superdirective afitenna 
having infinite gain. This undesirable result will be pre- 
vented by appropriately restricting the construction of the antenna. 

The antenna to be considered in this section consists of a rectangular hole in an infmite, 
perfectly conducting plane, which will be taken as the xy-plane. The origin of coordmates will 
be centered at the mouth of the hole, as shown in Fig. 23. The aperture dimensions, which are 
the length and width of the hole, are x„ wavelengths in the x-direction and y^ wavelengths in the 
y-direction. In some manner, which we will not investigate, electric fields w^ll be produced 
on the xy-plane at the mouth of the hole when a generator is connected to the antenna terminals. 
It will be assumed that all the power from this generator is radiated into the upper hemisphere. 

The gain of the antenna, the ratio of (15) to (16), is 

n21lP(i,m)l^|P(f.m)l" + |P,(f,ni)|"] 



n[lP^(f,m)P + lPy(f.m)|^ + lP,U.r«)l1 di dm 



The conducting plane imposes 
and E (x, y, o) are zero except for (x, y) 
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components P (i, m) and P (i, m), which are the Fourier transforms of the tangential fields, must 
therefore be well-behaved functions of direction. Although we would expect P^(i, m) to have a 
peak in the direction of maximum power density, the sharpness of that peak is limited by the 
aperture size. (It is easily shown, for example, that 


|P^(i,m)l^di dm 


is less than or equal to the aperture area expressed in square wavelengths.) T1 
gular plane-wave component, P (f, m), is linked to the other two by the divergei 
i P^ + m Py + n P^ = 0 . 


lird rectan- 
relation, 

(10a) 


Except possibly in the directions, n = 0, corresponding to the horizon, P^(i, m) will also be a 
well-behaved function of direction; its peak will again be limited in sharpness by the aperture 
dimensions. In other words, regardless of how the antenna is adjusted, its gain in any direction, 
except possibly toward the horizon, will remain finite. 

Because of the divergence relation (10a), only two of the three field components may be 
adjusted independently. It will be convenient to work with E^(x, y, o) and E^(x, y, o) or their 
transforms, P^(i, m) and P^d, m). The rectangular shape of the aperture has been chosen to 
simplify the calculations. Since either tangential field, E^(x, y, o), is zero outside the rectan- 
gular aperture, its transform, P,(i, m), according to the Sampling Theorem, is completely 
determined by a set of samples taken on a grid of points in the im-plane. Rows of sampling 
points are spaced by l/x^ in the i-direction and by l/y in the m-direction: 




and the Fourier series representing the corresponding tangential field is 

1 _ . , -j2Jr[x(i +r/x )+y(m +s/y )] 

5r^ I Pi<^o + ""o^ ®/yo> ® 

for - 1/2 x^ < X < 1/2 x^ and - 1/2 y^ < y < 1/2 y^ (92) 

0 otherwise . 

These relations may be derived in the same way as (59) and (60) were derived; in fact, they 
may be inferred from (59) and (60). The constants and m play the same role as the con- 
stants x^ and y^ did in (59) and (60); they may be chosen arbitrarily to displace the sampling 
grid with respect to the origin in the fm-plane (the zenith). It will often be convenient to pick 
the direction (f^, m^) to coincide with the direction of maximum power density. 

In maximizing the available signal power, we will adjust the coefficients, P (i + r/x 
“o * ^/*o' ™o^ Strictly speaking, there is a doubly infinite set of 

each but, as a practical matter, only a finite number of points (1 + r/x , m + s/y ) on the 
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sampling grid fall within f ^ + m^ = 1; Only a finite number of points correspond to real direc- 
tions. Changing the coefficients corresponding to samples taken very far outside of the unit 
circle will not have much effect on the antenna pattern in real directions, as can be seen from 
(91). Practically then, we may be content to adjust a finite number of coefficients. 

We will start with the expression (30) for the available power, which we will write in 
matrix notation as 


y dn{FX [«.] X F*]} 

^ ^ u.h. 

J dn{FX F*]} 
u.h. 

in which [*] = [$(i, m)] is the matrix of power densities (32), and 
Fjj(i, m) 

F] = F(l,m)] = 

Fy(f, m) 


The next step is to write the available power in terms of the x- and y-plane-wave components: 
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all quantities being functions of (i, m). If the matrices involved are extended, the integrand may 
be written in terms of the pattern samples, (91b), as 

p X 4N> X p*]- 

in which p]- is the column matrix of pattern samples, 

Pxll 


Pyl 

Py2 


NxxVl "xxV2 


Nxy¥l ^xy ¥2 ' 


Nyx¥l V'l'2- 


• • • 

V¥2 • 


The extra marks on the matrix brackets indicate the partitioning shown in (97) and (98). If the 
denominator integrand is extended in the same way, the available power may be written 
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j dn{p X -[Nh xp*]-} 


dn{p X xp*)-} 


In this form, only the square matrices -[N}- and -[D]- a 


^G^ = y ^N^ c 


the available power may finally be v 


It will be recalled from Sec. IV that the denominator of the expression for the available 
power is proportional to the radiation resistance of the antenna, and the numerator is propor- 
tional to the mean-square open-circuit voltage. The same is true of (100). This means that 
any matrix -[H]- of finite order is positive definite Hermitian, so that its inverse exists, and 
the matrix -[G}- of like order is at least non-negative definite Hermitian. [The denominator 
matrix -[H]- is actually symmetric, as may be seen from (95).] The expression (100) may be 
maximized by holding the denominator constant and maximizing the numerator. (This amounts 
to maximizing the mean-square open-circuit voltage while holding the radiation resistance 
constant.) If a procedure similar to the one leading to (74) is used, the condition for maximum 
available power is found to be 


^ X -[G]- = Q! ^ X -[H}- , (101) 

where a is again the Lagrangian multiplier. Post-multiplying both sides of (101) by p*]-, we 
have 

p X -[G]- X p*]- = a p X -[H). X p*)- , 

so that the maximum available power is just {A^/Zr,) a. The problem of solving (101) for its 
maximum characteristic value and corresponding vector (the set of pattern samples) is essen- 
tially identical with the general characteristic value problem (88) discussed in Sec. VI. The 
main difference is that we have not shown a constraint corresponding to (89). It will be noted 
that if pl- is multiplied by a constant in (100), that expression is not changed, so no constraint 
is necessary. On the other hand, if we had stated the available power to be 

^a = ^ P P*F 


the solution vector resulting from (101) would have to be normalized to satisfy the above c 
straint. These are just two equivalent ways of stating the problem. 
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t about the order of the G- and H-matrices (the 
idered. In general, if the pattern samples in- 
ints on the sample grid from which appreciable 
n available power cannot be increased appreciably 


)f the number of pattern sample 


The foregoing maximization is complicated by the fact that both the optimum pattern shape 
and the optimum polarization of the receiving antenna are determined simultaneously. As we 
discovered, it is really necessary to proceed in this manner if a general power-density matrix 
is under consideration. On the other hand, if the transmitting antenna is horizontally polarized, 
for example, it would appear to be obvious that the receiving antenna should also be horizontally 
polarized. We will next investigate the validity of such conjecture by seeking the conditions 
under which the optimum antenna, when driven at its terminals, produces tangential fields 
E^(x,y, o) and (x, y, o) that are in a fixed ratio, independent of the location of the point (x, y, o) 
on the aperture plane. Equivalently, we ask under what conditions the optimum antenna has the 
pattern P^(i, m) = a P^(f, m), with a some complex constant . 

If the power densities are sufficiently well-behaved, as they will be in practice, the pattern 


each direction (f, m), for if it dii 
In other words, if the optimum £ 


to the denominator integrand of (96) in 
£Ould be available if Py(f, m) ^ a P^(i, m). 
e pattern P (f , m) = a P (i, m), we must 


l,a X [N'f, m) ] = ^(f , m) ^ X [D(i, m) ] . 

We have written the condition under which (102) maximizes the ratio of the numerator integrand 
to the denominator integrand of (96). Note that the Lagrange multiplier p depends on the direc- 
tion (f, m). Post-multiplying by the transformation matrix [T], given by (38), and using (94) 

£uid (95), we have 


{ 1, a X [T]j } X[i] = [T]^" 

3 that the Lagrange multiplier is a characteristic 


characteristic value of the power-density matrix ii 
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direction of space. The corresponding characteristic vector (consisting of a pair of horizontal 
and vertical plane-wave components, incidentally) is proportional to 

a X [T] ^ ^ - n(m — ai), {1 + am) 

^/i^ + m^ 

and may be taken to be 

,n(m - ai), (f + am) 

We seek the conditions under which 


n(m - ai), (i + am) X 


with p an arbitrary (real, positive) function of (i, m) and a an arbitrary (complex) constant. If 
/3 is eliminated from the two equations, there results, after some algebra, 

(n^ |m-af|^-|i+ am|^) $jj.^ = n(m-a=f'f) (i+ am) , (104) 

the desired relation among the power densities that allows the polarization of the receiving 
antenna to be fixed in such a way that P^U. m) = a P^(f, m) before the pattern shape is determined. 

The condition (104) does not impose a particularly severe restriction on the power-density 
matrix. It is obviously satisfied by uniform noise, for which = 0 and $ = ■}> : This 

was, of course, to be expected, since the available noise power does not depend on the adjust- 
ment of the Mtenna. (104) is also satisfied by the power-density matrix 



[$(i,m)J = iif.m) 


n |m - af I n(m - a*f) (i + am) 

n(m - ai) (i + am*) |i + am|^ 


in which $(i, m) is an arbitrary real positive function of direction and a is an arbitrary complex 
constant. It will be noted that this matrix is singular, so that the polarization of the plane wave 
arriving from each direction (i. m) is also fixed in time. We have given examples of extreme 
cases in which (104) is satisfied. 

If (104) is satisfied, so that the antenna pattern is of the form (102), we have, after post- 
multiplying (103) by P*]. 

£ X [N] X P*] = /?(i, m) £ X [D] X P*] 

= ^(f.m) |P^(i.m)|2_^x [D]X M 

a*J 

= ^(i.m) |P^{i.m)|2 A(f,m) , 
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A(f, m) = [D] X J 

= (l-mV21mRe{a} + iai^l-,2) , (, 05 ^ 

o that the power available with optimum polarization, as expressed by (96), becomes 
j ^(^.m) |p^(i,m)|^ A(f,m) dn 

w = a! -6 

a 2r) 


m)|^ A{e, m) dB 


(106) 


with PU, m) the greater of the two characteristic values of [$(i, m)] in each direction of space. 

M tois point. It should be clear that we are justified in regarding m. m), or more properly 
(A /Zri) as Uie power density. 

Having determined the polarization of the antenna, it remains to maximize the available 
power by adjusting the shape of its pattern. We need work only with the x-samples, P (i +r/x , 

'"o'''®/yo'- are numbered consecutively as before to form the set 

Pi" 


and if we use the column matrix of corresponding f's given by (91c), 


then, from (91a), we can write 

|P (f,m)|^ = p xf] xj;,x P*] . 


] = y pU, m) f] X £ A(i, m) dfi 




As before, the condition for maximum available power is 

£ X [G] = a _£X [H] , (111) 

and that maximum power is just [A^/Zn] a. It may be noted that both integrands (109) are real, 
so that both matrices [G] and [H] are now real. This means that the characteristic vector p] 
will also be real; the conjugate signs may be removed from (110). 


A considerable simplification of the problem resulted when we were able to determine the 
optimum polarization of the antenna before finding its pattern shape; the order of the G- and 
H-matrices was cut in half. Further simplification is possible if the aperture dimensions of the 
antenna are large compared with the wavelength. In that case, the -functions given by 

(91c) vary rapidly with direction and tend to bunch in one direction. Let us examine the de- 
nominator integral (109), giving the element of [H] in the i* row and column. Let the orig- 
inal numbering (r, s) correspond to i and the original numbering (p, q) correspond to k; 

sin7r[x^(f-fj- r] sin 7r[y^(m - m^)- s] 

“ik - 7r[x^(f-i^)-r] 7r[y^(m - m^)- s] 


sinjr[x^(f-i^)-p] sinir[y^(m - m^)- q] 

7r[Xjj(f-i^)-p] 7r[yQ(m-m^)- q] 


dj dm 


For large x^ and y^, the behavior of the integrand is largely determined by the four _ 

functions, so that the remaining factor, A(f, m)/n. may be removed from the integral after being 
set equal to its value at the point of maximum contribution of the -functions. (The fact 

that l/n-*“ at the horizon does not cause any difficulty, since 

I -- 

u.h. u.h. 

exists.) In addition, the range of integration may be extended to cover the entire fm-plane. 

The double integral may then be split into the product of two single integrals of the form 



sinjr[x^(f-i^)- r] sinx[x^(f - f^)- p] 

’r[Xo(i-io)-r] >r[x^(f-fg)-p] 


\ P = 

0 otherwise 
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The denominator integral under consideration is therefore approximated by 
A(i^ + r/x^, m^ + s/y^) 


J 1 - (f^ + r/ x^)^ - (m^ + s/y^)^ 


if k = i and (f + r/x r + (m + s/y r < 1 


so that the matrix [H] of (109) becomes diagonal and of finite order. The same approximations, 
when applied to the more general denominator integral (99), result in a matrix of finite order 
having three non-zero diagonals; 


Cj 0 
0 


It is interesting to note that these same approximations can be deduced from a physical 
argument, which, although more difficult to justify than the preceding mathematical argument, 
provides a valuable insight into the meaning of the approximations. The denominator of the 
expression for available power, we recall, is proportional to the radiation resistance of the 
antenna, or, originally, the power emitted by the antenna when a unit current is applied to its 
terminals. The approach leading to the approximations consists of calculating the transmitted 
power from (17), using the Fourier series (92) for the fields. 

Strictly speaking, (92) applies only to the tangential fields, E^(x,y, o) and Ey(x,y, o), and 
only on the aperture plane of the antenna. On the other hand, any one of the terms of the sum 
behaves like a component of a plane wave over the aperture. It would seem reasonable that 
the normal field, E (x,y,o), should behave like the third component of a plane wave, the tan- 
gential components of which are given by (92), at least we would expect this behavior of the 
normal field at points within the aperture that are well removed from the edge. In other words, 
if the aperture is large enough, (92) is approximately valid for all three field components. 
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In order to mske use of (17), it is necessary to calculate 
ponents, which requires a knowledge of the field components 
vicinity of the aperture plane. Again it would seem reasonal 
term represented the component of a plane wave progressing 


:tend (92) as though each 
2 upper hemisphere. For 


^ 2 Pi<^o ■'/'"o' “o ®/yo> 


-j27r[x(f +r/x )+y(m^+s/y^)+s 


for (x, y) inside the aperture 


a = J 1 - (fjj + r/x^)^ - (m^ + s/yj 


If this approximation is used in calculating ( 


ismitted power becomes, after 


V = ^ ^ J 1 - (io + - (m^ + s/yy |P.(1^ + r/x^, m^ + s/y^)| 


Since the radical is pure imaginary outside of the unit circle, only a finite number of terms 
are involved in the double sum. If the divergence relation (10a) is used to express in terms 
of P and P , this expression for the transmitted power becomes identical with the expression 
that would be obtained with the approximate H-mati ix, (113): 

W = _P X X p*h . 

In other words, the physical argument leads to the approximation (113) for the H-matrix; it 
would also produce the diagonal matrix with coefficients given by (112) if the polarization of 
the antenna were predetermined. 

Another way of stating the result of the physical argument is that the transmitted power is 
calculated as though a collection of true plane waves were crossing the xy -plane and progress- 
ing into the upper hemisphere; the perfectly conducting sheet, however, is removed from the 
xy-plane, and only that power passing through the rectangular outline of the aperture is counted. 
In any case, it appears that the effects of the conducting edge of the aperture may be neglected 
if the aperture is large enough. Presumably, the edge effects of other large antennas can also 
be neglected. A dish, for example, which does not have an aperture bounded by the edge of a 
hole in an infinite conducting sheet, might be handled using the same approximations. Regard- 
less, of the structure of the antenna, its optimization will not make it superdirective if these 
approximations are used. 
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The same approximations, when applied to the numerator integrals (99) and (109), lead to 
numerator matrices -[GJ- and [G] of the same order as the corresponding denominator matrices, 
but the simplification of the integrals is not so spectacular in this case. The numerator ma- 
trices do not become diagonal because the power densities, as well as the -functions, vary 

rapidly witli direction, preventing the integrands from exhibiting the orthogonality possessed by 
the functions by themselves. 


Before closing this section, we will choose a representative power-density pattern and 
optimize the patterns of a number of antennas of different size for this power-density pattern. 

We are primarily interested in the shapes of the optimum patterns and in how much "gain loss" 
may be recovered by using the optimum antennas; we are not particularly interested in polariza- 
tion. Accordingly, we will suppose that (104) is satisfied, so that the polarization may be fixed 
in such a way that P^(f, m) = a P^(f, m). We will choose a power-density pattern centered on 
il , m ) and having the Gaussian shape 


Zri 


PU.m] = 




The parameters x^ and determine the sharpness of the power-density pattern; we will sup- 
pose that they are both so large that most of the power arrives from within a very small cone 
of directions centered on the real direction (f , m ). The power-density pattern has been nor- 
malized in such a way that the total power is 


u.h. 


di dm 


di dm = 1 watt 


As yet, we have not given a physical interpretation for this "total power," but one will be forth- 
coming. 

We will choose large antennas having aperture dimensions comparable to x^ and y^. The 
grid of points at which pattern samples are taken will be so fine and the number of pattern 
samples we will require will be so small that the factor A(f. m)/n can not only be removed 
from the matrix integrals.(109). but it can be set equal to its value in the direction of maxi- 
mum power density. Under these circumstances, there can be no doubt about the pattern 
samples all being taken in real directions. The denominator matrix, (109). will be 
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and the numerator matrix becomes 



f di dm 


The available power, (110), may be stated to be 
P, X [g] ^ Pi 


subject to the constraint 


with 


^ X p] = 1 

[g]=y]’ jLdidm . 


(116) 


(117) 


(118) 


The maximum available power is therefore just n^x^y^ where is the maximum charac- 
teristic value of [g]. 

Before calculating the elements of [g] for the intermediate cases, x^ x x^ and y^ » y^, we 
will examine the extremes. First, if the aperture is very large (x^ » x^ and y^ » y^), the 
antenna can resolve considerable detail in the power-density pattern; the maximum available 
power will be greatest in this case. We will suppose that the pattern samples are taken in 
directions reasonably close to The behavior of the integrand of (118) is then deter- 
mined primarily by the -functions, so that 



In this case, regardless of the antenna pattern [within reason: the pattern should still be con- 
centrated in the vicinity of (i^, m^)], the maximum available power will be n^x^y^/x. 

On the other hand, if the aperture is small, but still has dimensions of many wavelengths 
so that the approximations for the denominator matrix are valid, if 1 « x^ « x^ and 
1 « y « y^, the antenna cannot resolve the solid angle from within which most of the power 
density aurrives. As far as we could tell by adjusting the antenna pattern, the antenna is being 
illuminated by a point source. The integrand of (118) behaves like the power-density pattern. 

Throughout the present example we will use the renumbering scheme shown in Fig. 24 for 
the pattern samples, so that the first sample is taken in the direction of maximum power den- 
sity. The other elemental patterns, fj^, are all zero in that direction, as can be seen from 
(91c). In the present extreme case, then. 
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The integral is recognized as the "total power," 1// = 1 watt, (115), and the maximum available 
power is n x^y . The quantity n x y is the area in square wavelengths of the projection of the 
aperture onto a plane normal to the direction (i , m ). This makes the physical significance of 
the "total power" clear: It is the power per square wavelength of aperture area that would be 
available at the terminals of a large antenna if all the power density were concentrated in a 
single direction, and if that direction were normal to the aperture plane. 

The characteristic vector making the power n x y available from a point source is 


p] = 0 


which corresponds to the pattern, » 

sin IT x^(f - l^) sin ir y^(m - m^) J 

’rx^(f-fj xy^(m-m^) ' -j 

The field component that gives rise to this pattern is of constant amplitude and linearly varying 

phase over the aperture: ,| 

J -j2jr(xf^+ym^) 

x^y^ ® ■ -s 

E^(x, y, o) = ■ -1/2 x^ < X < 1/2 x^ and -l/2 y^ < y < l/2 y^ j 

ti 

0 otherwise . 
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If the pattern giving maximum plane-wave gain 
source is evidently n^x^y^^g^^, whereas the power 
total power would be i 


s used, the power available from the diffuse 
available from a point source of the same 


x°y . The "gain loss" is therefore the ratio of 1 


And since the 


maximum power available from the diffuse source 


3 n X y a .the ratio of g^^ to 


fraction of the "gain loss" that can be recovered by optimum pattern adjustment. 
It will be convenient to plot a normalized available power. 


a n^x^y^ 

a function of a normalized aperture dimension. 


When the aperture dimensions are small compared with x^ and y^, or if all the power density 
is concentrated in a single direction, this normalized available power may be obtained with the 
pattern giving maximum plane-wave gain and is equal to the normalized aperture area: w^ = 

(X /x.) (y /y,); when the aperture dimensions are of the order of x^ and y^, the normalized 
available “power is w- = (x^x^) (y„/y^) g^^ if the pattern giving maximum plane-wave gam is 
used, and it is = (x /x^) (y^y^) if the optimum pattern is used; when the aperture di- 
mensions are large compared with x, and y^, the normalized available power reaches the upper 
limit, w = l/ir. 

Next, we discuss the calculation of the elements of the g-matrix, (118). Let the original 
numbering (r, s) correspond to the new number i (Fig. 24), and let (p, q) correspond to k. The 
element of [g] in the i*’^ row and column is then 

CC ’*'1^1 "^*1 ^^1 ^ df dm 

Sik = J J “ " 

sinir[x^(i-fjj) -r] sinir[y^(m - m^) - s] sin ir[x^(i - - p] sin ir[y^(m - m^) - q] 

^ ir[x^(i-f^j) -r] ir[y^(m - m^^) - s] 7r[x^(i - t^) - p) ’r[yo(ni - m^) - q] 

The integral may be broken down into the product of two integrals of the same form: 

With x^(f — replaced by X, the first of these becomes 

with c the normalized aperture dimension (120). The evaluation of this integral is considered 
in the appendix. 
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With reference to (121), it should be noted that the new number i corresponds to the original 
numbering (r, s), not (r, p); it is not possible to express the matrix elements as the product 
^ik ~ ®’i\’ we discovered in Sec. VI, would lead to a particularly simple result. 

However, the fact that each matrix element can be expressed as a product of a factor depend- 
ing on x-dimensions and a factor depending on y-dimensions, (121), will simplify the calcula- 
tions considerably. (121) is a direct consequence of our having chosen a power-density pattern 
that can also be written as a product of two functions: In the narrow range of directions where 
(114) has appreciable magnitude, it may be written as the product of a function of the direction 
cosine t and a function of the direction cosine m. As a result of (121), the optimization of the 
antenna pattern may be split into two simpler one-dimensional problems, as we shall now see. 

The available power, (116), will be maximized subject to the constraint, (117), by the set 
of pattern samples satisfying 


t should be noted that the ri 


s not possible to express the 


^ik = ^i\’ a® ' 

However, the fact that ( 


ampling grid numbering. 


in terms of the original numbering. (The unfortunate double use 
cause any confusion.) Using (121), we have 



of the symbol p should not 


which may be rewritten in matrix notation as 

[Jit X [P] X [K] = n[p] , (123) 

in which [p] = [Pj.g] is the matrix of pattern samples, ar>'anged in a 

rectangular array corresponding to the array of points on the sampling grid. [J] and [K] are 
square, symmetric matrices, as may be seen from (122). If the sampling grid has R rows 
(different values of r) and S columns, [J] will be of order R and [K] will be of order S. 

Now suppose that the characteristic value problems, 

^X[J] = XU ] 

and V ^ (124) 

V X[K]= pv I 


ll 
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so that 

[P1=U]XV . . = 

is a solution to (123). In fact, there are R solutions to the first of equations (124) and S solutions 
to the second equation and hence RS solutions to (123) of the form (125). But, since there are 
RS points on the sampling grid and hence RS pattern samples to be determined, all solutions to 
(123) are of the form (125), with u_, v_. A, and h- determined from (124). 

If the constraints, 

^Xu] = l and ^ X v] = 1 (126) 

are satisfied by the characteristic vectors of (124), the constraint (117), which may be written 
variously as 

1 = £ X p] = ^ Pj^ = 2 Prs 
i r, s 

= 1 \ ^s = 2^1: '^s 
r, s r s 

is also satisfied. 

Physically, the characteristic vector £ determines the behavior of the antenna pattern as a 
function of the direction cosine i, and ^ determines its behavior with m; the optimum pattern 
may be expressed as the product of a function of i and a function of m, as may be seen from (91). 
The tangential fields over the aperture may be written in a corresponding product form, as may 
be seen from (92). The solutions of (124) corresponding to the greatest characteristic values, 

\ and ji, will result in the maximum available power. 

Since we have chosen a power -density pattern that is symmetric about 1 = 1^ (and about 
m = m ), we expect an optimum pattern that is also symmetric about t = (and about m = m^): 
We expect the elements of the characteristic vectors of (124) to satisfy u_^ = u^ (and v_^ = v^). 

If the sampling grid extends symmetrically about the point (r, s) = (0, 0), the elements of the 
characteristic vectors of (124) will actually show either even symmetry, u_^ = u^, or odd sym- 
metry, u = — u and u = 0, as we Will next demonstrate. 

We need the relations Jp^. = J^p and = Jj.p. which may be seen from (122). When these 

are used, the J-matrix becomes 
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The simplification ii 


• • (J22+J-22) <Jl2"''-12^ 


[M] X [J] X [M]^ = 


I-' 22 ^‘'-22' 


0 (Ji2“'^-12^ • 

° (■^1? " *'^22 “ ■^-22' • 


The original characteristic 
vector set ^ , given by 


istic value problem. (124), may be rewritten in terms of the transformed 


The transformed characteristic value problem reads 
^ X [M] X [J] = X [M] , 

or, after post-multiplication by [M]^, 

u; X{[M] X [J] X [M]^} = Xii: X [M] X [Ml^ 


since the transformation matrix is orthogonal, i.e. [M]^ = [M]''. The corresponding constraint. 
(126), reads 

^xu] = id X [M]x [M]^xu'] =u; xu'l = 1 • 


u 


The characteristic values of the simplified matrix, [M] X [J] X [M]^, are characteristic 
values of the sub-matrices located on its principal diagonal. The corresponding characteristic 
vectors, have zero elements in either the last (R - i)/Z positions or the first (R + l)/2 posi- 
tions. depending on whether the associated characteristic value is a characterisUc value of the 
upper left-hand sub-matrix or of the lower right-hand sub-matrix. If the characteristic value 
is from the upper left-hand sub-matrix, only the first (R + l)/2 rows of [M] are needed to calcu- 
late the characteristic vector u^ of the original problem (124), and it will be noted from [M] that 
in this case, u = u . On the other hand, if the characteristic value is from the lower right-hand 
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sub-matrix, only the last (R — l)/z rows of [M] are needed, and u = — u and u =0. Since 

Uq = 0 implies a null in the antenna pattern in the direction of maximum power density it is to 
be expected that the optimum pattern will show even symmetry, so that the upper left-hand 
sub-matrix will have the greatest characteristic value. This was found to be the case in all the 
calculations made. 

Although the discussion up to this point has presumed a rectangular aperture and a power- 
density pattern with unequal parameters, x^ and y^, the calculations were further simplified by 
making y^/y^ = = c, so that the normalized aperture dimensions are equal. It is conven- 

ient to regard the aperture as being a square with side and the power density pattern as being 
circularly symmetric about (i^,m^), so that y^ = x^. A square sampling grid with an odd num- 
ber of rows and columns has been used. These simplifications make [J] = [K], so that [p] is 
square; the maximum characteristic values and vectors of (124) are also equal; and 

At the start it was not clear how many samples would be necessary, so the greatest char- 
acteristic values of a 25 X 25 g-matrix and a 49 X 49 g-matrix were obtained. These corre- 
spond to sampling grids with 5 rows and columns and 7 rows and columns, or to J-matrices of 
order 5 and 7, respectively. In none of these calculations did the greatest characteristic value 
of the 25 X 25 g-matrix differ from the greatest 
characteristic value of the 49 X 49 g-matrix by 
more than 0.0047 decibels. We may easily con- 
clude that the sampling grid with 5 rows and 
columns is sufficient for the present example. 

The difference in characteristic vectors is 
just barely discernible in the patterns in direc- 
tions well removed from (i , m ). 

The normalized available power (119) is 
plotted against the normalized aperture di- 
mension (120) in Fig. 25. The bottom curve 
is the power available from an antenna having 
maximum plane-wave gain, and the sloped line 
is the power available from a point source of 
the same total power, y/ , when the same an- 
tenna is used. The decibel difference between j || 

these curves is the "gain loss." The middle ^ 

curve is the power available with optimum pat- j 

tern adjustment. Since a characteristic value : | 

problem must be solved to obtain each point on •: 

the curve, only the seven marked points have 
been obtained. The decibel difference between 
this curve and the lower curve is the recover- * 

able "gain loss." I' 
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Fig. 25. The available power increases with the 
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Fig. 26. The power advantage of an antenna having an optirr 
pattern is not very great . 
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Figures 27 through 30 are 
of the 49 X 49 matrices. The o 
ture dimension for which calcu) 
the power-density pattern has c 
the antenna patterns do not, bei 
cross sections of the patterns i 
the aperture dimension increas 
appear to remain about the san 
The -pattern giving maxii 


is disappointingly small; it is plotted in Fig. 26. 
set of antenna patterns plotted from the characteristic vectors 
imum patterns are for alternate values of the normalized aper- 
dons were made: c = l/l.5, l/0.7, l/0.3, 1/0.15. Although 
•cular symmetry about the direction of maximum power density, 
use the aperture is not circular. Figures 27 through 30 are 
the plane m = m . Also, the patterns really get narrower as 
5 , but plotting them against x^(f — l^) makes the antenna patterns 
width, while the power-density pattern appears to get broader, 
im plane-wave gain and the power-density pattern are plotted 


The required distribution of tangential elei 
excitation of the antenna at its terminals, is pi 
magnitude of electric field across the aperture 
over the aperture at rates proportional to the c 
plotted for the normalized aperture dimension, 
"gain loss" is the greatest. The normalizatioi 


1 of tangential electric field on the aperture plane, resulting from 
s terminals, is plotted in Fig. 31, which shows contours of constant 
TOSS the aperture. (In the present example, the phase varies linearly 
oportional to the direction cosines and m^.) Figure 31 has been 
erture dimension, c = l/0.7, the value at which the recoverable 
The normalization is such that 


yj' lE^(x,y,0)|^ 


We may conclude from the foregoing example that, if the power density pattern has a 
Gaussian shape, the effort required to design and construct an optimum antenna will generally 
not be sufficiently compensated by increased available power. In many cases, an antenna 
designed to have low side lobes will deliver more power than an antenna of the same size hav- 
ing maximum plane -wave gain. It seems reasonable to extrapolate these conclusions to any 
power -density pattern having a single, fairly well defined maximui 
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VIII. OPTIMUM DIVERSITY 



- 0.5 - 0.4 - 0.3 - 0.2 - 0.1 0 0.1 0.2 0.3 04 0.5 


Fig. 31 . Contours of electric field amplitude on the aperture plar 
for an antenna having an optimum pattern . 


In discussing the optimum coupling network for the set of electrically separate sources of 
Fig. 16, Sec. VI, we found that in the extreme case when the open-circuit source voltages are 
perfectly correlated, the coupling network can deliver a power to a single load equal to the sum 
of the powers available from each source separately. But if the source voltages are completely 
uncorrelated, the other extreme, the maximum power available from the coupling network, is 
only the greatest of the powers available from each source separately. We also remarked that, 
if the sources represent elements of a receiving array, the desirable perfect coherence between 
source voltages may be expected to persist for a period of time short compared with the fading 
period of the arriving signal, but that a time-varying coupling network would then be required 
to take advantage of this coherence. Although it will be convenient, initially, to regard the time- 
varying coupling network as an extension of the work of Sec. VI, the resulting apparatus is really 
a form of spaced-antenna diversity which results in the maximum signal-to-noise (power) ratio 
of the resultant signal. 

The signal-to-noise ratio available from a number of coherent but noisy sources has been 
maximized by Kahn, this maximization has been established in a more general form by 
Brennan, and implementations of the mathematical result have been described by Kahn and by 
Mack, who calls his equipment "combiner diversity." But whereas the equipment discussed 
by Kahn and Mack combines the signals after detection, we will arrange to combine signals before 
detection; The equipment suggested by this section will continue to deliver the maximum signal- 
to-noise ratio of the resultant signal even if the signal power available from each antenna is less 
than the noise power. 

We will again represent the terminal-pair behavior of an array of N antennas by the set of 
disconnected sources on the left in Fig. 16. This representation entails no loss in generality, 
for, if there is appreciable mutual impedance between the terminal pairs of an actual array, the 
array elements may be interconnected as in Figs. 20 and 21 to produce an equivalent set of N dis- 
connected sources, capable of supplying the same power. The analysis of this section will then 
apply to the equivalent set of sources. 

In terms of the matrix of transmission coefficients, t], and the matrix of normalized open- 
circuit signal voltages of the sources in Fig. 16, V/’/SR] (Eq. 75), the instantaneous signal power 


only the greatest of the powei 


By "instantaneous" power we mean the power averaged over an integer number of radio-frequency 
cycles, but not averaged for such a long time that the modulating waveform and particularly the 
fading waveform has a chance to change. As we discovered in Sec. VI, this instantaneous power 
is maximized subject to the constraint, 

_;^Xt*] = l , (68a) 

by the set of transmission coefficients satisfying 
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h e an arbitrary phase shift, and the maximum instantaneous power li 


(79a) 


which is the sum of the instantaneous signal powers available from each source separately. If 
the transmission coefficients are continuously adjusted according to (80a) to maintain the maxi- 
mum instantaneous power (T9a), the mean signal power is again the sum of the mean signal powers 
available from each source separately. The mean noise power, as discussed in Sec. VI, is inde- 
pendent of the adjustments of the transmission coefficients if, as we shall assume, the noise 
arrives uniformly from all directions in space. 

It is enlightening to regard this continuously adjusted array as a single composite antenna. 

We note that the maximum power available from an array (of electrically separate elements) 
illuminated by a point source is also the sum of the powers available from each element. Thus, 
if the indi\T.dual elements of our composite antenna are not subject to "gain loss" when they are 
iUuminated by a diffuse source, the composite antenna itself will not suffer "gam loss" when 
illuminated by the same diffuse source'. 

To witMnthe approximations discussed in the last section, which are valid when aperture 
dimensions are large compared with the wavelength, the power available from an element of the 
composite a.n-enr.a is proportional to its aperture area if the elemental antenna is too small to 
■resolve' the power-density pattern, but approaches an upper limit as the aperture area is in- 
creased (see Fig. 25). Apparently, if its total aperture area is fixed, the composite antenna will 
deliver the greatest power if its elemental antennas are reduced in size (and increased in number) 
until they can no longer "resolve* the power-density pattern. From the point of view taken in the 
introductory section, a scatter link employing this continuously adjusted array will be most eco- 
nomical of the apermre areas of transmitting and receiving antennas if the transmitting antenna 
and each of the elemental antennas of the receiving array are not able to resolve the scattering 
region. 

We have seen that a coctinuously adjusted array can deliver more signal power to a matched 
load than a single large antenna having the same aperture area. But if each transmission coef- 
ficient of the time-varying coupling network is to be made proportional to the conjugate of the 
signal \-oltage at the terminals of the corresponding array element, as required by (80a), it would 
appear that the signal at each array element would have to be known exactly in order for the com- 
posite signal to be received with maximum signal-to-noise ratio! In addition, we should be con- 
cerned about whether the time variations of the transmission coefficients will perturb the com- 
posite signal. Fortunately, for one choice of the "arbitrary phase shift," G of (80a), a complete 
knowledge of the signals is not required, and the time variations of the transmission coefficients 
do not perturb the composite signal, as we shall now see. 

It will be necessary to differentiate between signal variations caused by modulation of the 
transmitter and signad variations caused by fading. As shown in Fig. 32. we will suppose that 
the transmitted signal is modified differently in following each of the paths from transmitting 
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antenna to the various receiving array elements. M is the vector transmitted signal having 

e lan carrier requency and modulation described by the complex time function M. (The 
instantaneous current at the transmitting antenna terminals would be proportional to Re {U 
In traversing the path to the 

in proportion to the complex time function K.. which may be regarded as the gain of the path. 
Before the signal appears as the voltage at the terminals of the i*'’ array element, noise is added; 
the various path noises will be assumed to be uncorrelated and equal in power. 

Note that the model of Fig. 32 precludes the 
effects of multipath transmission: If the trans- 
mitter is pulsed, a single replica of that pulse 
is received by each of the elemental antennas. 
Although this assumption is more restrictive 
than the assumptions discussed in Sec. II, it is 
often met in practice and is necessary for the 
proper functioning of the apparatus we will de- 
vise. Since the path gains are assumed to 
be independent of frequency (over the band of 





ysRi 


v^^Rn 


Fig. 32, Model of the radio paths from the transmitting luuepenuem oi irequency (over the bai 
cel^n^orU'" """ '-'I— — P-d by the transmitter) 


provision is made for the time-of-flight of the 
o - of no consequence here. 

From Fig. 32, the matrix of normalized open-circuit signal voltages of the array elements 
may be written 


signal from transmitter to receiving site; this 


xTsr 


(127) 


in which fC is the matrix of path gains. The maximum instantaneous power, (79a), becomes 

W^. = iMl^KXK*] , ( 128 ) 

and it may be obtained if the transmission coefficients, from (80a), are 


M* 






(Ml. J0\ ^ 


(129) 


The corresponding normalized open-circuit signal voltage at the output of the coupling network 
is then 


V \ 

— ( = I /-y- in the notation of Fig. 16; 

\ / 

= t = |m| ej^'/KXK*] , (130) 

From (130), we see that the signal emerging from the coupling network is proportional, 
among other factors, to the root-sum-square path gain - a real, slowly varying function of time. 
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* * nha<!e of the output signal would also be con- 

If we chose the phase shift 0 to be a cons an , ^ proportion to the magnitude of the 

stant; only the amplitude of the output signa w . coupling network as a carr 

modulation. A frequency-modulated signal would emerge from th^e coup 
fading Slowly in amplitude onl,.. ^^der the same cond - , 

--- 

choice is 

ie _ M 

e = Arg M ; e - ■ 

Por this choice, the transmission coefficients need only be adiusted in accordance with the si 
varying path gains: 

,K*, (12' 

n/KXK*T 

the output signal is directly proportional to the modulation: 


and neither of the difficulties appears. 

If the transmitter is unmodulated (M = 1), the instantaneous power available from 

varying coupling network, (1Z8), is 


KXK*1 = E lK.|^ . 

i=l 

unless the magnitudes of the path gains are perfectly correlated, the relative spread in the dis- 
tribution of power output due to fading will be less (in a mean-square sense) than the relative 
spread in the distribution of power available from one of the array elements. In other words, 
the decibel fading range is reduced by the continuously adjusted array - another of its advantages. 
Staras^* discusses the case in which the jK, 1 are independent and Rayleigh-distributed. These 
conditions would be expected to hold for short periods of time, say of the order of an hour, and 
if the array elements were appreciably separated in space. He points out that the power output 
follows a Chi-square distribution (having 2 N degrees of freedom), which he plots for N up to 10. 

It remains for us to find an appropriate 

I K, NOISE FRONT IDENTICAL mcasure of the path gains and to show how these 

1 - ^^ 1 filters , ^ 

1 ^ I 1 ^ ^ may be obtained and used in an electronic reaii- 

r i . zation of the time-varying coupling network. 

; 1 I ‘ ■ First, since the instrumentation is to be elec- 

1 JL I 1 - 

© * ~T •'09 ' ' tronic, we will need workable signal levels. 

: : : : and it will be expedient to filter out the noise 

I rQ •-© 1 I —— except within the band of frequencies occupied 

by the transmitter. Suitable arrangements are 
Fig. 3l The model has been modified slightly to in- indicated in the block diagram of Fig. 33, which 
elude appropriate amplifiers and filters. 
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will be used as a model of the problem at hand. The outputs of this model, the E's, are the am- 
plified and filtered array element voltages. It will be assumed that the front-end gains have been 
set so that the noise levels of these outputs are the same, and, as discussed previously, the 
noises in the different outputs are uncorrelated. These considerations are important, as we shall 
see. If the front ends include mixers, as will generally be convenient, it is also important that a 
common local oscillator be used. The signal-to-noise ratios of the outputs will all be degraded 
by the noise figures of the front ends, but this unavoidable loss would also appear if a passive 
coupling network and single front end were used. 

After amplification, the question of obtaining maximum power does not arise; it is the signal- 
to-noise ratio that must be maximized. The latter problem, however, is not new; it has already 
been solved: We found that, in addition to providing maximum power in a matched load, the cou- 
pling network output. 



also results in a maximum signal-to-noise ratio when the transmission coefficients are adjusted 
according to (129a), since the output noise power is independent of the adjustment of the trans- 
mission coefficients (when the noises from the separate sources are uncorrelated and equal in 
power). By setting the front-end gains in Fig. 33 to equalize the output noise levels, we have 
brought the output signal voltages, the E's, into proportion with the normalized open-circuit 
signal voltages of the array elements, i. e.. 



Thus, the signal-to-noise ratio of the combined output, 

Eo = £,XE1 , (132) 

or a multiple thereof, will also be a maximum. 

It will be convenient to lump the fixed, equal front-end gains, 

Ei 

vT ’ 

with the path gains, K.. in Fig. 33, so that the output signal voltages may be written 

E=MK . (133) 

Defined in this way, the path gains of Fig. 33 are no longer equal to, but are in direct proportion 
to, the path gains of Fig. 32. This will make no difference in the "transmission coefficients," as 
may be seen from (129a), and, of course, the absolute level of the combined output (132) is im- 
material: We will continue to use the same symbols, K^. 
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The combined output. 


K*XE] 



then, would have the maximum signal-to-noise ratio, and it wo uld have a 

root-mean-square signal level would be proportional to ^/KX K*], as may be 


(132a) 

fixed noise level. The 
seen from (13 3) and (132a). 


The combined output, jj 

Fq = k;!:xE] = 2 K?Ej 


(134) 


would also have the m aximum signal-to-noise ratio, but the root-mean-square noise level would 
vary in proportion to n/K X K*], and the root-mean-square signal level would be proportional to 
K X K*]. (The "mean" in these statements is interpreted as an appropriate short-time average.) 
to the design of a receiver, it is often convenient to use automatic volume control to hold the 
root-mean-square signal level fixed. This may be accomplished after the outputs of Fig. 33 
have been combined. It will be simplest to realize a circuit having the combined output (134), 
for which the combining coefficients are simply the conjugates of the path gains. 

We return to the problem of finding an appropriate measure of the path gains, which are 
defined by (133). If it were not for the noise, we might measure M at the transmitter and E^ at 
the receiver, giving K. = E./M. This possibility is not too far-fetched, for we hope to recover 
M at the receiver as well. It is the noise that causes the real difficulty; this estimate, at any 
instant of time, would be intolerably noisy. We recall that the path gains, which are the signals 
received when the transmitter is unmodulated (M = 1), vary only slowly; they describe the fading. 
Thus a short-time average, taken, at each instant, over a period of the past that is small com- 
pared with the fading period, would provide an accurate estimate of a path gain and at the same 
time would reduce the noise power of an instantaneous estimate by a factor equal to the product 
of the averaging time and the bandwidth of one of the filters shown in Fig. 33. For example, if 
the system (and fUter) bandwidth were 100 kcps and the fading were slow enough so that we could 
average over 0.1 second of the past, the averaging process would reduce the noise power by a 
factor of 10,000. Using this averaging process, we could obtain a continuous good estimate of 
a path gain for output signal-to-noise ratios of that path greater than about -40 decibels . This 
short-time average will be indicated by the notation } - 

If the short-time average is performed by a linear passive network, as we intend, contribu- 
tions to the average will be weighted in accordance with how far back in the past they occurred. 
Another kind of weighting may be considered: A reasonable measure of a path gain would be 



with the arbitrary time function, g. real and non-negative. To be consistant, the same weighting 
fxmction, g, would have to be used in measuring the gain of each path. 

It would be sensible to pick g in such a way that it increases with the magnitude of the mod- 
ulation: g a monotonic increasing function of |M |. This kind of choice gives more weight to the 
instantaneous estimate, E^/M, when the signal amplitude, and hence the signal-to-noise ratio. 


no 
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:~2Z T:~- TT " 

t • J u .. mpared with the maximization of signal-to-noise ratio ob- 
tained when the outputs of Fig 33 are u- 

First the / • „ • P-'operly combined, (134), using the measured path gains, 

r 1 ^ --y estimate. K.. but neglecting 

sincTthl rat T indistinguishable from the fading, 

the shirt performing 

the hort-time average. We may therefore disregard the variations in the estimated path gains 

due to noise, and we are left with the prospect of having a reduced signal-to-noise ratio result- 
ing from improperly determined, but momentarily fixed "transmission coefficients." Except 
when the estimates are violently in error, this effect is again small. As we discovered in Sec. VI 
the correct estimates result in a stationary point of the signal-to-noise ratio. Thus, a small frac- 


tional error. 



in one of the estimates would result in a smaller fractional error of the order of | e a reduc- 
tion of course, in the signal-to-noise ratio. 

A reasonable choice of weighting function is g = | M |. We will accept this choice not because 
It IS optimum in some sense, but because it leads to particularly simple circuitry. Our estimate 


Not having a copy of the transmitted signal at the receiver, we will use the phase of the com- 
bined output (134) as an estimate of the phase of the transmitted signal. This estimate will be 
quite satisfactory as long as the signal-to-noise ratio of the combined output is only slightly 
greater than unity, since errors in this estimate affect only the measured path gains. Thus, our 
measurement of path gain will give 



and we will obtain the combined output 


N 


JArgF^ 


I 


(134a) 


Containing, as it does, the output F^^ on both sides, (134a) is reminiscent of the equation of 
an oscillator. In fact, the circuit for obtaining this output, shown symbolically in Fig, 34, is an 
oscillator of sorts . The inputs to this circuit are the outputs of the front ends and filters shown 
in Fig. 33; the output appears at the upper right of Fig. 34. The tuned circuits marked per- 
form the short-time averages. They are all tuned to some (radian) frequency and have 

bandwidths of the order of the reciprocal of the averaging time (10 cps in the numerical example 
previously mentioned). The limiter delivers a constant- amplitude output sinusoid (chosen to be 



UNCLASSIFIED 



one volt in Fig. 34). and the phase of its output is the same as the phase of the input. The limiter 
should have sufficient dynamic range to limit properly on noise. 

Noting that the instantaneous frequencies 
of the signals emerging from the tuned circuits 
must be close to coj, we see that the instantaneous 
frequency of the signal leaving the second mul- 
tiplier in each channel must be about — oi^. 

(It is assumed that product components at the 
sum frequency are suppressed.) The output, 
which is the sum of the outputs of the second 
multipliers, may then be represented as having 
the carrier frequency as has been 

done in Fig. 34. It may be seen, by following 
an assumed output through the processes of 
limiting, first multiplication, averaging, sec- 
ond multiplication, and summing, that the cir- 
cuit of Fig. 34 does indeed produce the output 
given by (134a). 

The properties of this combining circuit 
may be deduced either from the figure or from 
(134a). Suppose, first, that the phases of all 



Fig. 34. Anelectronicreaiization of tF 
coupling network. 


! time-varying 


the inputs are advanced suddenly by the same amount, simulating a step increase in the phase 
of the modulation. Nothing will happen immediately to the tuned-circuit voltages (the short- 
time averages will remain fixed), but the phases of all the second-multiplier outputs and hence 
the phase of the output will advance by the same amount. The phases of the first-multiplier 
outputs (the quantities being averaged), however, will not change. Hence there will be no addi- 
tional, long-term change in the output phase. In other words, the phase of the modulation is 
reproduced in the output, perturbed, possibly, by noise and fading. The output behaves differ- 
ently with amplitude changes, however. For rapid changes in the amplitude of the modulation, 
the output follows the input because the tuned circuits cannot follow the input changes. But if 
the modulation amplitude changes slowly, at a rate less than or comparable to the fading rate, 
the output aunplitude varies as the square of the modulation amplitude: The tuned circuits can 
follow these slow changes, and the magnitude of the loop gain of this "oscillator" is twice mul- 
tiplied by the amplitude of the input signals. From another point of view, the decibel fading 
rauige in the outjmt is twice the range that would be obtained from a circuit using normalized 
transmission coefficients and producing the combined output (132a). As we have already noted, 
this difficulty may be mitigated by an automatic volume control operating on the combined out- 
put Fq. 

The behavior of the circuit of Fig. 34 in the absence of signal and under the influence of 
noise alone has not been studied, either experimentally or analytically. Nevertheless, we may 
conjecture about its behavior when the inputs are filtered Gaussian noise. First, if the noise 
inputs are uncorrelated and equal in power, we would expect the tuned- circuit outputs to vary 
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slowly in amplitude and phase, behaving very nearly like independent filtered Gaussian noises 
having the same power. The application of even a very weak signal should cause the tuned- 
circuit outputs to tend to fall in a fixed proportion in such a way as to maximize the signal-to- 
noise ratio at the output. 

On the other hand, if one of the noise powers were greater than the others, we would expect 
the corresponding tuned-circuit output also to be large. In fact, we would expect the large tuned- 
circuit output to suppress the others in such a way that the combined output follows the greatest 
input noise rather closely. Similarly, if the input noises were partially correlated, we would 
expect the tuned-circuit outputs to tend to adjust themselves so as to maximize the output noise. 

In either of these conditions, the combining circuit may be regarded as biased; its adjustment 
favors one set of input signals. A fairly strong signal would then be required to overcome the 
bias and cause the tuned-circuit outputs to readjust themselves properly for that signal. 

The circuit of Fig. 34 bears a resemblance to the circuits devised by Price^^ for another 
application. The resemblance would be stronger if the limited output were replaced by a stored- 
reference signal, having some of the characteristics of the signal expected from the transmitter. 
If a stored-reference signal is available at the receiver, such an arrangement is preferable to 
the circuit of Fig. 34 when the output signal-to-noise ratio is small. 
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■^rp = I du y n(v- |) n (V + |) dv 



Fig. A-2 . The two (l-fijncfions are symmetri- 
cally displaced. 


j_ p{c-|ul 

^%^i(c-|u| 


We will call the second integral 1^. .-Apparently 
it is an even function of u; the magnitude of its 
integrand is the product of the two n-functions 
shown in Fig. A-2. Note that if |u| > c, = 0, 
so we need only be concerned with the range 

■ 

and if p r, 

. N, 


■(c-|u|)/2 


" l“l (r-p) _ ^ (r - pf s 


• |u| (r-p) 


TT I u| (r - p) 


since r and p are integers. 

Taking the special case p = r first. 




L 
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If this integral is split into the sum of two integrals having tl 
and if u is replaced by -u in the first, the sum may be writt« 


J^^.2 (1-x) e-<^^">'cos 2.r 


In the more general case when p r, 


-(ffu)^ ^|u| (r-p) p-j7ru(r+p)/c 

FITTF^ ® 


The range of this integral may also be cut in half as befoi 


2(-l)t^-P> r -(^v 

IT c ( r - p) 


u (r - p) TT u (r + p) j 


and. with the aid of the appropriate trigonometric identity, 


rp TT (r - p) r p' 


■ C sin 27r u ^ du 


The lnt.gr.1. (A-i) end (A-3) c.nnot b. performed by el.me»..r, methods. However rf 
the upper limit, are .rbltrarll, estend.d to they cm. both be e.alu.t.d In term, of tabulated 
functions, and If the parameter c (the normalized aperture dlm.n.l.nl is large, the error 
,ol.ed will be small. Even with c a. small as unity, the eaponcn.ial has dropped to about 
» V .0-‘ at the present upper limit'. Fortunately, we wlU be interested m value, of c of he 
order of unity or greater, so the approalm.tlon appears to be reasonable. We wdl calculate 
approximation, for the integrals and give an upper bound for the errors Involved. St.rtmg with 
the second integral, (A-3), 

, . sin 1 . rx dx . dim I f eJ^'" dxl . 
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Completing the square, 
I “ I 


"{r 


g-[(7rcx)^-j2irrx-(r/c) ] ^ 
^-(?rcx-jr/c)^ 


“} ■ 


)r, with the complex variable, z = (jt cx- j r/c). 


|-3i-3574| 


Fig . A-3 . Path of integration i n the cc 
plex plane . 




.-(r/c)' 


Ir “ ^ e"'^''"' <£m 

2 pr/c 2 


0 follow the imagina 
A-3. The integral < 
aginary axis. With 


The path of integration may be taken to follow the imaginary axis to the origin and then along the 
real axis to infinity, as shown in Fig. A-3. The integral along the real axis is pure real, so we 
are left with the integral along the imaginary axis. With z = -jy, the latter integral becomes 


The function 


D(x) = e‘* j* e^ dy 

s tabulated. In terms of this function, the approximation to reads 
Ir ” ^ D(|) . 

The approximation to is 

-incx)^ ^^3 27 r rx dx 


(1 - x) e 
„-(>rcx)^ 




-(JTCX) 


COS Zir rx dx 


^13 + 14 • 

The first integral may be written 

I 3 = dx , 

and it may be performed by completing the square. The result is 
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integral may be integrated by parts; 


- C cos Zir rx d [e 


The integral is recognized as the approxim 


la = ^ ■ 


and the approximation to 


The approximation to J^p, from (A-2) and (A-5), 


rp~ ^V-P)c " 

;t may be noted from (A-1) that the central integral, can be expressed exactly 11 




This central integral is of primary importance in determining the power available at the antenna 
terminals; the power available from an antenna with a square aperture having maximum plane- 
wave gain is proportional to g,, = The exact value (A-S) is used in the numerical calcu- 
lations described in Sec. VII. „ 

we will not attempt to relate the error in the maximum characteristic value of the g-matnx 
to the errors involved in the approximations (A- 6 ) and (A-7). In giving an upper bound for the 
latter, it will be appropriate to compare the errors with the central integ 
1 j^ (exact) - J^.p (approx.) | 
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and it will be convenient to use 


°° n/7c 

for c greater than 1/2, this error in the error calculation will reduce the tru< 
lying between 0.6 and 1.0. 

From (A-1), the fractional error in J is 


c ^ (x- 1) e 


-(ttcx) 


xdx 




2x dx- 2V7 c e-<^=*) 



Fig. A-4. The fractional errors soon become insignificant as the nor- 
malized aperture dimension increases. 
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error by a factor 


<A-9) 


--i 
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And from (A-2) and (A-3), the fractional error ii 






Zit px] dx ] 




(A-IO) 


The upper bounds on the fractional errors, (A-9) and (A-IO), are plotted in Fig. A-4. It will be 
noted that the fractional errors are less than 1 part in 10^ if c is greater than 0.66. The ap- 
proximations (A-6) and (A- 7) are only used in cases where c >0.66. 
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